IDEMPOTENT/TROPICAL ANALYSIS, THE HAMILTON-JACOBI AND 

BELLMAN EQUATIONS 

GRIGORY L. LITVINOV 

^Nj , Abstract. Tropical and idempotent analysis with their relations to the Hamilton-Jacobi and matrix 

Bellman equations are discussed. Some dequantization procedures are important in tropical and idem- 

C " 3 ' potent mathematics. In particular, the Hamilton-Jacobi-Bellman equation is treated as a result of the 

^\J ' Maslov dequantization applied to the Schrodinger equation. This leads to a linearity of the Hamilton- 

Jacobi-Bellman equation over tropical algebras. The correspondence principle and the superposition 
j^ , principle of idempotent mathematics are formulated and examined. The matrix Bellman equation and 

its applications to optimization problems on graphs are discussed. Universal algorithms for numerical 
algorithms in idempotent mathematics arc investigated. In particular, an idempotent version of interval 

— ^ ^ analysis is briefly discussed. 
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In dear memory of my beloved wife Irina. 



1. Introduction 

In these lecture notes we shall discuss some important problems of trop- 
ical and idempotent mathematics and especially those of idempotent and 

^ ■ tropical analysis. Relations to the Hamilton-Jacobi and matrix Bellman 
equations will be examined. Applications of general principles of idempo- 

P ■ tent mathematics to numerical algorithms and their computer implemen- 
tations will be discussed. 

Tropical mathematics can be treated as a result of a dequantization of 
the traditional mathematics as the Planck constant tends to zero taking 
imaginary values. This kind of dequantization is known as the Maslov de- 
quantization and it leads to a mathematics over tropical algebras like the 
max-plus algebra. The so-called idempotent dequantization is a generaliza- 
tion of the Maslov dequantization. The idempotent dequantization leads 
to mathematics over idempotent semirings (exact definitions see below in 
sections 2 and 3). For example, the field of real or complex numbers can 
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be treated as a quantum object whereas idempotent semirings can be ex- 
amined as "classical" or "semiclassical" objects (a semiring is called idem- 
potent if the semiring addition is idempotent, i.e. x © x = x), see [32H12]. 
Some other dequantization procedures lead to interesting applications, e.g., 
to convex geometry, see below and |iHl[5Hl [5B]. 

Tropical algebras are idempotent semirings (and semifields). Thus trop- 
ical mathematics is a part of idempotent mathematics. Tropical algebraic 
geometry can be regarded as a result of the Maslov dequantization ap- 
plied to the traditional algebraic geometry (O. Viro, G. Mikhalkin), see, 
e.g., P21[721[751IMH^ . There are interesting relations and applications to 
the traditional convex geometry. 

In the spirit of N.Bohr's correspondence principle there is a (heuristic) 
correspondence between important, useful, and interesting constructions 
and results over fields and similar constructions and results over idempotent 
semirings. A systematic application of this correspondence principle leads 
to a variety of theoretical and applied results J32H13] , see Figure [H 

The history of the subject is discussed, e.g., in [32], with extensive bib- 
liography. See also |I5l El IHl 1201 1221 HDHlZiaS]. 



V.P. Maslov's idempotent superposition principle means that many non- 
linear problems related to extremal problems are linear over suitable idem- 
potent semirings. The principle is very important for applications including 
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Figure 1. Relations between idempotent and traditional mathematics. 
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numerical and parallel computations. See V.P. Maslov's original formula- 
tion in P3HS5], as well as P[Ill[I5l[I7l[IHl[2ni[22lE3lEM3lll5], and below. 



2. The Maslov dequantization 

Let R and C be the fields of real and complex numbers. The so-called 
max-plus algebra Rmax = RU{— oojis defined by the operations x (By = 
max{a:, y} and x Qy = x + y. 

The max-plus algebra can be seen as a result of the Maslov dequantization 
of the semifield R+ of all nonnegative numbers with the usual arithmetics. 
The change of variables 

X \-^ u = hlogx, 

where h > 0, defines a map ^h- R+ — > R U {— oo}, see Fig. [21 Let the 
addition and multiplication operations be mapped from R+ to RU {— oo} 
by $/i, i.e. let 

u(BhV = hlog{exp{u/h) + exp{v/h)), uQ v = u + v, 

= -oo = $/,(0), 1 = = $/,(!). 

It can be easily checked that w ©/^ -u — > max{w, -u} as /i — > 0. This 
deformation of the algebraic structure borrowed from R+ brings us to the 
semifield Rmax, known as the max-plus algebra, with zero = — oo and 
unit 1 = 0. 

The semifield Rmax is a typical example of an idempotent semiring; this 
is a semiring with idempotent addition, i.e., x(Bx = x for arbitrary element 
X of this semiring. 

The semifield Rmax is also called a tropical algebra.The semifield R*^'^^ = 
$/i(R+) with operations (Bh and (i.e.+) is called a subtropical algebra. 

The semifield Rmin = R U {+00} with operations © = min and = + 
(0 = +00, 1 = 0) is isomorphic to Rmax- 

The analogy with quantization is obvious; the parameter h plays the role 
of the Planck constant. The map x 1— > |x| and the Maslov dequantization 
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for R+ give us a natural transition from the field C (or R) to the max-plus 
algebra Rmax- M/^e will also call this transition the Maslov dequantization. 
In fact the Maslov dequantization corresponds to the usual Schrodinger 
dequantization but for imaginary values of the Planck constant (see below). 
The transition from numerical fields to the max-plus algebra Rmax (or 
similar semifields) in mathematical constructions and results generates the 
so called tropical mathematics. The so-called idempotent dequantization is 
a generalization of the Maslov dequantization; this is the transition from 
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Figure 2. Deformation of R+ to R*^'''. Inset: the same for a small value of h. 
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basic fields to idempotent semirings in mathematical constructions and 
results without any deformation. The idempotent dequantization generates 
the so-called idempotent mathematics, i.e. mathematics over idempotent 
semifields and semirings. 

Remark. The term 'tropical' appeared in [SD] for a discrete version of 
the max-plus algebra (as a suggestion of Christian Choffrut). On the other 
hand V.P. Maslov used this term in 80s in his talks and works on economical 
applications of his idempotent analysis (related to colonial politics). For 
the most part of modern authors, 'tropical' means 'over Rmax (or Rmin)' 
and tropical algebras are Rmax and Rmin- The terms 'max-plus', 'max- 
algebra' and 'min-plus' are often used in the same sense. 

3. Semirings and semifields. The idempotent correspondence 

principle 

Consider a set S equipped with two algebraic operations: addition © and 
multiplication 0. It is a semiring if the following conditions are satisfied: 

• the addition © and the multiplication are associative; 

• the addition © is commutative; 

• the multiplication is distributive with respect to the addition ©: 

and 

{xey)Qz = {x(Dz)e{yQz) 
for all x,y, z G S. 

A unity (we suppose that it exists) of a semiring S is an element 1 E S 
such that iQx = xQl = X for all X E S. A zero (if it exists) of a semiring 
S is an element E S such that 7^ 1 and 0©a: = x, O0a: = a:0O = Ofor 
all X E S. A semiring S is called an idempotent semiring if. x (B x = x for 
all X E S. A semiring S with neutral element 1 is called a semifield if every 
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nonzero element of S is invertible with respect to the multiphcation. For 
the theory of semirings and semifields the reader is referred, e.g., to |I26]. 

The analogy with quantum physics discussed in Section E] and below 
leads to the following idempotent correspondence principle: 

There is a (heuristic) correspondence between important, useful and in- 
teresting constructions and results over the field of complex (or real) num- 
bers (or the semifield of nonnegative numbers) and similar constructions 
and results over idempotent semirings in the spirit of N. Bohr's correspon- 
dence principle in quantum theory j^U\^^. 



This principle can be also applied to algorithms and their software and 
hardware implementations. Examples are discussed below; see also |59] - li^ 

4. Idempotent analysis 

Idempotent analysis deals with functions taking their values in an idem- 
potent semiring and the corresponding function spaces. Idempotent anal- 
ysis was initially constructed by V. P. Maslov and his collaborators and 
then developed by many authors. The subject is presented in the book of 
V. N. Kolokoltsov and V. P. Maslov ||33] (a version of this book in Russian 
was published in 1994). 

Let S be an arbitrary semiring with idempotent addition © (which is 
always assumed to be commutative), multiplication 0, and unit 1. The 
set S is equipped with the standard partial order :<: by definition, a ^ b 
if and only if a © 6 = 6. If S contains a zero element 0, then all elements 
of S are nonnegative: ^ a for all a E S. Due to the existence of this 
order, idempotent analysis is closely related to the lattice theory, theory of 
vector lattices, and theory of ordered spaces. Moreover, this partial order 
allows to model a number of basic "topological" concepts and results of 
idempotent analysis on the purely algebraic level; this line of reasoning 
was examined systematically in [|39j- [|57j and J18j . 
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Calculus deals mainly with functions whose values are numbers. The 
idempotent analog of a numerical function is a map X —)■ S^ where X is 
an arbitrary set and S is an idempotent semiring. Functions with values 
in S can be added, multiplied by each other, and multiplied by elements 
of S pointwise. 

The idempotent analog of a linear functional space is a set of S'-valued 
functions that is closed under addition of functions and multiplication of 
functions by elements of 5, or an ^-semimodule. Consider, e.g., the S- 
semimodule B{X, S) of all functions X — > 5 that are bounded in the sense 
of the standard order on S. 

If S = Rmax, then the idempotent analog of integration is defined by the 
formula 



1) I{(p) = / (p{x) dx = sup (p{x] 

Jx xex 



where (p G B{X^S). Indeed, a Riemann sum of the form '^<-p{xi) ■ di 

i 

corresponds to the expression 0(/?(a:i) (d oi = max{(y5(xi) + 0"^}, which 

tends to the right-hand side of ([ID as cr^ —> 0. Of course, this is a purely 
heuristic argument. 

Formula ([I]) defines the idempotent (or Maslov) integral not only for 
functions taking values in Rmax, but also in the general case when any of 
bounded (from above) subsets of S has the least upper bound. 

An idempotent (or Maslov) measure on X is defined by the formula 
m^(y) = supi/j{x), where i/j G B{X,S) is a fixed function. The integral 

x&Y 

with respect to this measure is defined by the formula 



(2) I^{(p) = / (p{x)dm^= / (p{x) Q tlj{x) dx = sup {(p{x) Qtlj{x)). 
Jx Jx xex 
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Obviously, if 5 = Rmm, then the standard order is opposite to the con- 
ventional order <, so in this case equation (E]) takes the form 

r® r® 

/ ip{x) dm^ = / (p{x) Q ip{x) dx = m[{(p{x) Qip{x)), 
Jx Jx ^ex 

where inf is understood in the sense of the conventional order <. 

We shall see that in idempotent analysis measures and generalized func- 
tions (versions of distributions in the sense of L. Schwartz) are generated 
by usual functions. For example the ^-functional 6y : (f{-) i— > (f{y) is 
generated by the function 

1^0, itx^y. 
It is clear that 

if{y) = / 6y{x) (p{x)dx = sup{6y{x) (p{x)) 

Jx X 



5. The superposition principle and linear equations 

5.1. Heuristics. Basic equations of quantum theory are linear; this is 
the superposition principle in quantum mechanics. The Hamilton-Jacobi 
equation, the basic equation of classical mechanics, is nonlinear in the con- 
ventional sense. However, it is linear over the semirings Rmax and Rmin- 
Similarly, different versions of the Bellman equation, the basic equation 
of optimization theory, are linear over suitable idempotent semirings; this 
is V. P. Maslov's idempotent superposition principle, see [63-65]. More 
generally, the idempotent superposition principle means that although 
some important problems and equations (related to extremal problems, 
e.g., optimization problems, the Bellman equation and its instances, the 
Hamilton-Jacobi equation) are nonlinear in the usual sense, they can be 
treated as linear over appropriate idempotent semirings. For instance, 
the finite-dimensional stationary Bellman equation can be written in the 
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form X = i7 X F, where X, i7, F are matrices with coefficients in 
an idempotent semiring S and the unknown matrix X is determined by 
H and F, see below and P [HllSllSni 1221 [2Sl 122]. In particular, standard 
problems of dynamic programming and the well-known shortest path prob- 
lem correspond to the cases S = Rmax and S = Rmin, respectively. It is 
known that principal optimization algorithms for finite graphs correspond 
to standard methods for solving systems of linear equations of this type 
(i.e., over semirings). Specifically, Bellman's shortest path algorithm cor- 
responds to a version of Jacobi's algorithm. Ford's algorithm corresponds 
to the Gauss-Seidel iterative scheme, etc. |I11II5]- 



The linearity of the Hamilton-Jacobi equation over Rmin and Rmax, which 
is the result of the Maslov dequantization of the Schrodinger equation, is 
closely related to the (conventional) linearity of the Schrodinger equation 
and can be deduced from this linearity. Thus, it is possible to borrow 
standard ideas and methods of linear analysis and apply them to a new 

Consider a classical dynamical system specified by the Hamiltonian 

N 2 

where x = (xi, . . . , xn) are generalized coordinates, p = [pi, . . . ^p^) are 
generalized momenta, rrii are generalized masses, and V{x) is the potential. 
In this case the Lagrangian L{x,x^t) has the form 

N .n 



L{x, ±, t) = ^ m^^ - V{x) 



i=l 

where x = (ii, . . . , ijv), Xi = dxi/dt. The value function S{x,t) of the 
action functional has the form 

S= L{x{t),x{t),t)dt, 

Jto 
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where the integration is performed along the actual trajectory of the sys- 
tem. The classical equations of motion are derived as the stationarity 
conditions for the action functional (the Hamilton principle, or the least 
action principle). 

For fixed values of t and ^o and arbitrary trajectories x{t), the action 
functional S = S{x{t)) can be considered as a function taking the set of 
curves (trajectories) to the set of real numbers which can be treated as 
elements of Rmm- In this case the minimum of the action functional can be 
viewed as the Maslov integral of this function over the set of trajectories or 
an idempotent analog of the Euclidean version of the Feynman path inte- 
gral. The minimum of the action functional corresponds to the maximum 
of e"*^, i.e. idempotent integral /r„at/is| e~^'^^''^^^D{x{t)} with respect to the 
max-plus algebra Rmax- Thus the least action principle can be considered 
as an idempotent version of the well-known Feynman approach to quantum 
mechanics. The representation of a solution to the Schrodinger equation 
in terms of the Feynman integral corresponds to the Lax-Olemik solution 
formula for the Hamilton-Jacobi equation. 

Since dS/dxi = pi, dS/dt = —H{p,x), the following Hamilton-Jacobi 
equation holds: 

Quantization leads to the Schrodinger equation 

(4) --^ = ^^ = //(^,x,)^, 

I ot 

where ?/' = ip{x^t) is the wave function, i.e., a time-dependent element 
of the Hilbert space L^(R^), and H is the energy operator obtained by 
substitution of the momentum operators Pi = j-£r and the coordinate op- 
erators Xi'. ip \-^ Xiip for the variables pi and xi in the Hamiltonian function, 
respectively. This equation is linear in the conventional sense (the quan- 
tum superposition principle). The standard procedure of limit transition 
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from the Schrodinger equation to the Hamilton-Jacobi equation is to use 
the following ansatz for the wave function: 'ijj{x^t) = a(a:, t)e*'^^^'^^/^, and 
to keep only the leading order as /i — > (the 'semiclassical' limit). 

Instead of doing this, we switch to imaginary values of the Planck con- 
stant h by the substitution h = ih, assuming h > 0. Then the Schrodinger 
equation (jll) becomes similar to the heat equation: 



(5) h^^H^-hg-,i.ju, 

where the real- valued function u corresponds to the wave function ^JJ. A 
similar idea (a switch to imaginary time) is used in the Euclidean quantum 
field theory; let us remember that time and energy are dual quantities. 

Linearity of equation ([Ij) implies linearity of equation (E]). Thus if ui and 
U2 are solutions of (E]), then so is their linear combination 

(6) u = XiUi + X2U2. 

Let S = hlnu or u = e^/^ as in Section [2] above. It can easily be checked 
that equation ([5]) thus turns to 

ds .,_ ^1 fdSV , A I d^S 






(7) _ 11-/ 

Thus we have a transition from (E]) to (Ej) by means of the change of vari- 
ables il) = e^l^ . Note that \'\\)\ = e^^^l^ , where Re^S is the real part of S. 
Now let us consider ^ as a real variable. The equation ([71) is nonlinear in 
the conventional sense. However, if S\ and 5*2 are its solutions, then so is 
the function 

(8) 5 = Ai ^iS/jAa ^2 

obtained from (0) by means of the substitution S = hXnu. Here the 
generalized multiplication coincides with the ordinary addition and the 
generalized addition ©/i is the image of the conventional addition under 
the above change of variables. As /i — > 0, we obtain the operations of the 
idempotent semiring Rnmx, i-e., © = max and = +, and equation ([7j) 
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becomes the Hamilton- Jacobi equation (jSj), since the third term in the 
right-hand side of equation (I7j) vanishes. 

Thus it is natural to consider the limit function 5 = Ai 5*1 A2 <5'2 
as a solution of the Hamilton-Jacobi equation and to expect that this 
equation can be treated as linear over Rmax- This argument (clearly, a 
heuristic one) can be extended to equations of a more general form. For 
a rigorous treatment of (semiring) linearity for these equations see, e.g., 
[331 US, [H5] . Notice that if h is changed to —h, then we have that the 
resulting Hamilton-Jacobi equation is linear over Rmin- 

The idempotent superposition principle indicates that there exist impor- 
tant nonlinear (in the traditional sense) problems that are linear over idem- 
potent semirings. The idempotent linear functional analysis (see below) 
is a natural tool for investigation of those nonlinear infinite-dimensional 
problems that possess this property. 

5.2. The Cauchy problem for the Hamilton-Jacobi equations. A 

rigorous "idempotent" appproach to the investigation of the Hamilton- 
Jacobi equation was developed by V.N. Kolokoltsov and V.P. Maslov j32] (a 
Russian version of this book was pubhshed in 1994); see also [TTIIHBII^I^ . 
Let us consider, inspired by a long tradition, the well-known Cauchy 
problem for the Hamilton-Jacobi equation (^. Given the action function 
at time T 

(9) ^(r, x) = St{x) = (p{x), X e R^, 

the Cauchy problem asks to reconstruct S{t, x) for x G R^ during the time 
interval <t<T. 

We shall discuss the min-plus linearity of this problem and denote by 
Ut the resolving operator, i.e. the map which assigns to each given St{x) 
the solution S{t, x) of the Cauchy problem in the interval < t < T. 
Then the map Ut, for each t, is a linear (over Rmin) operator in the space 
LSC(R'^, Rniin) of lower semicontinuous functions taking their values in 
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Rniin- Moreover Ut is an integral operator (in the sense of idempotent 
mathematics) of the form: 

(10) {Ut^){x)= ip{y)Kt{x,y)dy = mf{ip{y) + Kt{x,y)}, 

J y 

where Kt{x, ?/), as a function of ?/ G R", is bounded from below and lower 
semicontinuous. See [321IS5] for details. 

The operator Ut (as well as other integral operators, see Section [7] below) 
has the following property: 

(11) Ut{^if,) = ^{Utcp,), 

V V 

where {}Pv] is a bounded set of elements in LSC(R",Rniin)- So if we 
have such a family of functions SyiT^x^ and S{T,x) = / Si,{T^x)dv = 
infi,(5'i,(r, x)), then the solution of the Cauchy problem is expressed as 

S{t,x) = mf.jy{Sjy{t^X)). 

Relations between the "idempotent approach", viscosity solutions and 
minimax solutions in the sense of A.I. Subbotin JSSIISS] are examined, e.g., 
E^ in details; see also W.M. McEneaney [71]. To this end, let us mention 



m 



that more general Hamiltonians of the form H = H{t,x,p) (satisfying 
some additional conditions) and different kinds of solution spaces are also 
considered in the literature. 

The situation is similar for the Cauchy problem for the homogeneous 
Hamilton- Jacobi equation 

where H : H^ i— )► R is a convex (not strictly) first order homogeneous 
function 

H{p)= sup (f-p + g), /gR", geR, 
and y is a compact set in R'^+^. See [33]. 
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To develop a rigorous "idempotent" approach to differential equations 
and other problems, one needs an idempotent version of analysis and, es- 
pecially, functional analysis. See Section [7] below. 

6. Convolution and the Fourier-Legendre transform 

Let G be a group. Then the space B{G, Rmax) of all bounded functions 
G — > Rmax (see above) is an idempotent semiring with respect to the 
following analog ® of the usual convolution: 

{(p{x) ®tl)){g) == / (p{x) Q tlj{x^^ ■ g) dx = sup{(p{x) + tl){x^^ ■ g)). 

Jg xeG 

Of course, it is possible to consider other "function spaces" (and other 

basic semirings instead of Rmax)- 

Let G = R", where R" is considered as a topological group with respect 

to the vector addition. The conventional Fourier-Laplace transform is 

defined as 

(12) ip{x) ^ ip{0 = f e'^-'^ifix) dx 

Jg 

where e*^'^ is a character of the group G, i.e., a solution of the following 
functional equation: 

f{x + y) = f{x)f{y). 
The idempotent analog of this equation is 

fix + y) = fix) fiv) = fix) + fiy), 

so "continuous idempotent characters" are linear functionals of the form 

X \-^ ^ ■ X = ^iXi + • • • + ^nXn- As a rcsult, the transform in ( TT^ assumes 

the form 

/•© 

(13) (fix) \-^ (fi^) = / ^ • X Q (fix) dx = suipi^ • X -^ (fix)). 

Jg xeG 

The transform in ( TT3I ) is the Legendre transform (up to some change of no- 
tation) [53] ; transforms of this kind establish the correspondence between 
the Lagrangian and the Hamiltonian formulations of classical mechanics. 
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The Legendre transform generates an idempotent version of harmonic anal- 
ysis for the space of convex functions, see, e.g., [ jHT] . 

Of course, this construction can be generalized to different classes of 
groups and semirings. Transformations of this type convert the generalized 
convolution ® to the pointwise (generalized) multiplication and possess 
analogs of some important properties of the usual Fourier transform. 

The examples discussed in this sections can be treated as fragments of 
an idempotent version of the representation theory, see, e.g., [SD]. In par- 



ticular, "idempotent" representations of groups can be examined as rep- 
resentations of the corresponding convolution semirings (i.e. idempotent 
group semirings) in semimodules. 

7. Idempotent functional analysis 

Many other idempotent analogs may be given, in particular, for basic 
constructions and theorems of functional analysis. Idempotent functional 
analysis is an abstract version of idempotent analysis. For the sake of 
simplicity take S = Rmax and let X be an arbitrary set. The idempotent 
integration can be defined by the formula (1), see above. The functional 
/((/?) is linear over S and its values correspond to limiting values of the 
corresponding analogs of Lebesgue (or Riemann) sums. An idempotent 
scalar product of functions (p and i/j is defined by the formula 

{(fjt/j) = / (p{x) Q tlj{x) dx = sup {(p{x) Q tlj{x)). 
J X xex 

So it is natural to construct idempotent analogs of integral operators in 

the form 

/•© 

(14) ip{y) ^ {K(p){x) = K{x,y)Qip{y)dy = sup{K{x,y) + ip{y)}, 

Jy yeY 

where (p{y) is an element of a space of functions defined on a set y, and 
K{x,y) is an ^-valued function on A x y. Of course, expressions of this 
type are standard in optimization problems. 
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Recall that the definitions and constructions described above can be ex- 
tended to the case of idempotent semirings which are conditionally com- 
plete in the sense of the standard order. Using the Maslov integration, one 
can construct various function spaces as well as idempotent versions of the 
theory of generalized functions (distributions). For some concrete idem- 
potent function spaces it was proved that every 'good' linear operator (in 
the idempotent sense) can be presented in the form (Ej); this is an idem- 
potent version of the kernel theorem of L. Schwartz; results of this type 
were proved by V. N. Kolokoltsov, P. S. Dudnikov and S. N. Samborskii, 
I. Singer, M. A. Shubin and others. So every 'good' linear functional can 
be presented in the form (/? i— )■ {(p^ilj)^ where (, ) is an idempotent scalar 
product. 

In the framework of idempotent functional analysis results of this type 
can be proved in a very general situation. In [ J17] - I5UI[5^[S7] an algebraic 
version of the idempotent functional analysis is developed; this means 
that basic (topological) notions and results are simulated in purely al- 
gebraic terms (see below). The treatment covers the subject from basic 
concepts and results (e.g., idempotent analogs of the well-known theo- 
rems of Hahn-Banach, Riesz, and Riesz-Fisher) to idempotent analogs of 
A. Grothendieck's concepts and results on topological tensor products, nu- 
clear spaces and operators. Abstract idempotent versions of the kernel 
theorem are formulated. Note that the transition from the usual theory 
to idempotent functional analysis may be very nontrivial; for example, 
there are many non-isomorphic idempotent Hilbert spaces. Important re- 
sults on idempotent functional analysis (duality and separation theorems) 
were obtained by G. Cohen, S. Gaubert, and J. -P. Quadrat. Idempo- 
tent functional analysis has received much attention in the last years, see, 
e.g., PIISlEBHSniESlEB], 133]- W\ and works cited in |3S]. AU the results 
presented in this section are proved in [IS] (Subsections 17.11 - 17.41 ) and 
in [53 (Subsections O- inn!) 
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7.1. Idempotent semimodules and idempotent linear spaces. An 

additive semigroup S with commutative addition © is called an idempotent 
semigroup if the relation x © x = a: is fulfilled for all elements x E S. If 
S contains a neutral element, this element is denoted by the symbol 0. 
Any idempotent semigroup is a partially ordered set with respect to the 
following standard order: x ^ y if and only if x(By = y. It is obvious that 
this order is well defined and x (B y = sup{x, ?/}. Thus, any idempotent 
semigroup is an upper semilattice; moreover, the concepts of idempotent 
semigroup and upper semilattice coincide, see |[ID]. An idempotent semi- 
group S is called a- complete (or algebraically complete) if it is complete 
as an ordered set, i.e., if any subset X in S has the least upper bound 
sup(X) denoted by ©X and the greatest lower bound inf(X) denoted by 
AX. This semigroup is called b-complete (or boundedly complete), if any 
bounded above subset X of this semigroup (including the empty subset) 
has the least upper bound ©X (in this case, any nonempty subset Y in 
S has the greatest lower bound AY and 5 in a lattice). Note that any 
a-complete or 6-complete idempotent semiring has the zero element that 
coincides with ©0, where is the empty set. Certainly, a-completeness 
implies the 6-completeness. Completion by means of cuts |[ID] yields an 
embedding 5 —> 5 of an arbitrary idempotent semigroup S into an a- 
complete idempotent semigroup S (which is called a normal completion of 

S); in addition, S = S. The 6-completion procedure S ^ Sb is defined 
similarly: if 5 3 oo = sup 5, then Sb =S; otherwise, S = SbU {oo}. An 
arbitrary 6-complete idempotent semigroup S also may differ from S only 
by the element oo = sup S. 

Let S and T be 6-complete idempotent semigroups. Then, a homomor- 
phism / : 5 — > T is said to be a b-homomorphism if /(©X) = ©/(X) 
for any bounded subset X in S. If the 6-homomorphism / is extended to 
a homomorphism S —> T of the correesponding normal completions and 
/(©X) = ©/(X) for all X C S, then / is said to be an a-homomrphism. 
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An idempotent semigroup S equipped with a topology such that the set 
{s G iSls ^ 6} is closed in this topology for any 6 G 5 is called a topological 
idempotent semigroup S. 

Proposition 1. Let S be an a-complete topological idempotent semigroup 
and T be a b-complete topological idempotent semigroup such that, for any 
nonempty subsemigroup X in T, the element ©X is contained in the topo- 
logical closure of X in T. Then, a homomorphism f : T ^ S that maps 
zero into zero is an a-homomorphism if and only if the mapping f is lower 
semicontinuous in the sense that the set {t G T\f{t) :< s} is closed in T 
for any s E S. 

An idempotent semiring K is called a-complete (respectively b-complete) 
if K is an a-complete (respectively 6-complete) idempotent semigroup and, 
for any subset (respectively, for any bounded subset) X in i^ and any 
k e K, the generalized distributive laws /c0(©X) = ^{kQX) and (©X)0 
k = ©(X /c) are fulfilled. Generalized distributivity implies that any 
a-complete or 6-complete idempotent semiring has a zero element that 
coincides with ©0, where is the empty set. 

The set R(max, +) of real numbers equipped with the idempotent ad- 
dition © = max and multiplication = + is an idempotent semiring; in 
this case, 1 = 0. Adding the element = — oo to this semiring, we obtain 
a 6-complete semiring Rmax = R- U {— oo} with the same operations and 
the zero element. Adding the element +oo to Rmax and asumming that 
O0(+oo) = and a:0(+oo) = +oo for x 7^ and a: ©(+00) = +00 for any 
x, we obtain the a-complete idempotent semiring Rmax = R-max U {+00}. 
The standard order on R(max, +), Rmax and Rmax coincides with the or- 
dinary order. The semirings R(max, +) and Rmax are semifields. On the 
contrary, an a-complete semiring that does not coincide with {0, 1} cannot 
be a semifield. An important class of examples is related to (topological) 
vector lattices (see, for example, |ID] and pB], Chapter 5). Defining the 
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sum x©?/ as supjx, y} and the multiplication as the addition of vectors, 
we can interpret the vector lattices as idempotent semifields. Adding the 
zero element to a complete vector lattice (in the sense of IIDIES]), we 
obtain a 6-complete semifield. If, in addition, we add the infinite element, 
we obtain an a-complete idempotent semiring (which, as an ordered set, 
coincides with the normal completion of the original lattice). 

Important definitions. Let V be an idempotent semigroup and K 
be an idempotent semiring. Suppose that a multiplication k^x \-^ k Q x 
of all elements from K by the elements from V is defined; moreover, this 
multiplication is associative and distributive with respect to the addition 
in V and 1 Q x = x^ Q x = for all x E V. In this case, the semigroup 
V is called an idempotent semimodule (or simply, a semimodule) over K. 
The element Ov^ G V is called the zero of the semimodule V if kQOy = Oy 
and Ov^ © X = a: for any k E K and x E V. Let F be a semimodule over a 
6-complete idempotent semiring K. This semimodule is called b- complete 
if it is 6-complete as an idempotent semiring and, for any bounded subsets 
Q in i^ and X in V, the generalized distributive laws {(BQ)Qx = (B{QQx) 
and k (©A) = ©(/c A) are fulfilled for all k E K and x E X. This 
semimodule is called a-complete if it is 6-complete and contains the element 
oo = sup V. 

A semimodule V over a 6-complete semifield K is said to be an idempotent 
a-space {b-space) if this semimodule is a-complete (respectively, 6-complete) 
and the equality (AQ) 0a: = A{Q Q x) holds for any nonempty subset Q 
in K and any x E V, x y>^ oo = sup V. The normal completion y of a 6- 
space V (as an idempotent semigroup) has the structure of an idempotent 
a-space (and may differ from V only by the element oo = supV). 

Let V and W be idempotent semimodules over an idempotent semiring 
K. A mapping p : V ^ W is said to be linear (over K) if 

p{x © ?/) = p{x) © p{y) and p{k Q x) = k Q p{x) 
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for any x^y G V and k E K. Let the semimodules V and W be 6- 
complete. A linear mapping p : V ^ W is said to be b-linear if it is 
a 6-homomorphism of the idempotent semigroup; this mapping is said to 
be a-linear if it can be extended to an a-homomorphism of the normal 
completions V and W. Proposition 7.1 (see above) shows that a-linearity 
simulates (semi) continuity for linear mappings. The normal completion K 
of the semifield i^ is a semimodule over K. li W = K, then the linear 
mapping p is called a linear functional. 

Linear, a-linear and 6-linear mappings are also called linear, a-linear and 
b-linear operators respectively. 

Examples of idempotent semimodules and spaces that are the most im- 
portant for analysis are either subsemimodules of topological vector lat- 
tices [HS] (or coincide with them) or are dual to them, i.e., consist of linear 
functional subject to some regularity condition, for example, consist of 
a-linear functionals. Concrete examples of idempotent semimodules and 
spaces of functions (including spaces of bounded, continuous, semicontin- 
uous, convex, concave and Lipschitz functions) see in [331 SHI HE 157] ^^^ 
below. 

7.2. Basic results. Let V be an idempotent 6-space over a 6-complete 
semifield K, x E V. Denote by x* the functional V ^ K defined by the 
formula x*{y) = A{k G K\y :< kQx}, where y is an arbitrary fixed element 
from V. 

Theorem 2. For any x E V the functional x* is a-linear. Any nonzero 
a-linear functional f onV is given by f = x* for a unique suitable element 
xeV. IfK^ {0,1}, then X = e{y E V\f{y) ^ 1}. 

Note that results of this type obtained earlier concerning the structure of 
linear functionals cannot be carried over to subspaces and subsemimodules. 

A subsemigroup VF in V closed with respect to the multiplication by 
an arbitrary element from K is called a b-subspace in V if the imbedding 
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W —> V can be extended to a 6-linear mapping. The following result 
is obtained from Theorem [2] and is the idempotent version of the Hahn- 
Banach theorem. 

Theorem 3. Any a-linear functional defined on a b-subspace WinV can 
be extended to an a-linear functional on V. If x^y G V and x ^ y, then 
there exists an a-linear functional f on V that separates the elements x 
andy, i.e., f{x) ^ f{y). 

The following statements are easily derived from the definitions and can 
be regarded as the analogs of the well-known results of the traditional 
functional analysis (the Banach-Steinhaus and the closed-graph theorems). 

Proposition 4. Suppose that P is a family of a-linear mappings of an 
a-space V into an a-space W and the mapping p : V ^ W is the pointwise 
sum of the mappings of this family, i.e., p{x) = sup{pa{x)\pa G P}- Then 
the mapping p is a-linear. 

Proposition 5. Let V and W be a-spaces. A linear mapping p : V ^ W 
is a-linear if and only if its graph T in V x W is closed with respect to 
passing to sums (i.e., to least upper bounds) of its arbitrary subsets. 



In IIBI the basic results were generalized for the case of semimodules over 



the so-called reflexive 6-complete semirings. 

7.3. Idempotent 6-semialgebras. Let i^ be a 6-complete semifield and 
A be an idempotent 6-space over K equipped with the structure of a semir- 
ing compatible with the multiplication K xA ^ Aso that the associativity 
of the multiplication is preserved. In this case, A is called an idempotent 
b-semialgebra over K. 

Proposition 6. For any invertible element x E A from the b-semialgebra 
A and any element y G A, the equality x*{y) = !*(?/ x^^) holds, where 
l^A. 
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The mapping AxA^K defined by ttie formula (x, y) i— > (x, y) = l*(a:0 
y) is called the canonical scalar product (or simply scalar product). The 
basic properties of the scalar product are easily derived from Proposition E] 
(in particular, the scalar product is commutative if the 6-semialgebra A 
is commutative). The following theorem is an idempotent version of the 
Riesz-Fisher theorem. 

Theorem 7. Let a b-semialgebra A be a semifield. Then any nonzero a- 
linear functional f on A can be represented as f{y) = {y^x), where x G A, 
X ^ and (•, •) is the canonical scalar product on A. 

Remark 8. Using the completion precedures, one can extend all the results 
obtained to the case of incomplete semirings, spaces, and semimodules, 
see 



Example 9. Let B{X) be a set of all bounded functions with values belong- 
ing to R(max, +) on an arbitrary set X and let B{X) = B{X) U {0}. The 
pointwise idempotent addition of functions (ipi © ^2){x) = ^i{x) (f2{x) 
and the multiplication [cpi <-P2){x) = {(pi{x)) {(p2{x)) define on B{X) 
the structure of a 6-semialgebra over the 6-complete semifield Rmax- In 
this case, 'i-*{<-p) = sup^.^^^ (yi)(a:) and the scalar product is expressed in 
terms of idempotent integration: {(pi^(p2) = ^'^Pxexi'^ii^) © '^2{x)) = 

^'^Pxexi'^ii^) + '^2{x)) = J{(pi{x) <-P2{x)) dx. Scalar products of this 

X 

type were systematically used in idempotent snslysis. Using Theorems [2] 
and [7| one can easily describe a-linear functionals on idempotent spaces in 
terms of idempotent measures and integrals. 

Example 10. Let X be a linear space in the traditional sense. The idem- 
potent semiring (and linear space over R(max, +)) of convex functions 
Conv(X, R) is 6-complete but it is not a 6-semialgebra over the semifield 
K = R(max, +). 
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Any nonzero a-linear functional / on Conv(A, R) has the form 

(p \-^ f{(p) = sup {(p{x) + iIj{x)} = / ip{x) Qi/j{x) dx, 

X J X 

where ?/' is a concave function, i.e., an element of the idempotent space 
Conc(A, R) = - Conv(A, R). 

7.4. Linear operator, 6-seniiniodules and subsemimodules. In what 
follows, we suppose that all semigroups, semirings, semifields, semimod- 
ules, and spaces are idempotent unless otherwise specified. We fix a basic 
semiring K and examine semimodules and subsemimodules over K . We 
suppose that every linear functional takes it values in the basic semiring. 
Let V and W be 6-complete semimodules over a 6-complete semiring K . 
Denote by Li{V, W) the set of all 6-linear mappings from V to W . It is 
easy to check that Lb(y, W) is an idempotent semigroup with respect to 
the pointwise addition of operators; the composition (product) of 6-linear 
operators is also a 6-linear operator, and therefore the set L^iV, V) is an 
idempotent semiring with respect to these operations, see, e.g., [49]. The 
following proposition can be treated as a version of the Banach-Steinhaus 
theorem in idempotent analysis (as well as Proposition S] above). 

Proposition 11. Assume that S is a subset in Lfj{V, W) and the set {g{v) \ 
g G S} is bounded in W for every element v G V; thus the element f{v) 
= supg^g g{v) exists, because the semimodule W is b-complete. Then the 
mapping v i— > f{v) is a b-linear operator, i.e., an element of Lfj{V, W). The 
subset S is bounded; moreover, sup^S = /. 

Corollary 12. The set Lb{V, W) is a b-complete idempotent semigroup 
with respect to the (idempotent) pointwise addition of operators. If V = 
W, then Li,{V, V) is a b-complete idempotent semiring with respect to the 
operations of pointwise addition and composition of operators. 

Corollary 13. A subset S is bounded in Li,{V, W) if and only if the set 
{g{'^) \ g ^ S} is bounded in the semimodule W for every element v eV . 
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A subset of an idempotent semimodule is called a subsemimodule if it is 
closed under addition and multiplication by scalar coefficients. A subsemi- 
module y of a 6-complete semimodule W is b-dosed if V is closed under 
sums of any subsets of V that are bounded in W. A subsemimodule of 
a 6-complete semimodule is called a b-subsemimodule if the corresponding 
embedding is a 6-homomorphism. It is easy to see that each 6-closed sub- 
semimodule is a 6-subsemimodule, but the converse is not true. The main 
feature of 6-subsemimodules is that restrictions of 6-linear operators and 
functionals to these semimodules are 6-linear. 

The following definitions are very important for our purposes. Assume 
that W is an idempotent 6-complete semimodule over a 6-complete idem- 
potent semiring K and V is a subset of W such that V is closed under 
multiplication by scalar coefficients and is an upper semilattice with re- 
spect to the order induced from W. Let us define an addition operation 
in V by the formula x (B y = sup{a:,?/}, where sup means the least upper 
bound in V. If il' is a semifield, then y is a semimodule over K with 
respect to this addition. 

For an arbitrary 6-complete semiring K^ we will say that V is a quasisub- 
semimodule of VF if F is a semimodule with respect to this addition (this 
means that the corresponding distribution laws hold). 

Recall that the simbol A means the greatest lower bound (see Subsection 
7.1 above). A quasisubsemimodule V of an idempotent 6-complete semi- 
module W is called a /\- subsemimodule if it contains and is closed under 
the operations of taking infima (greatest lower bounds) in VF. It is easy to 
check that each A -sub semimodule is a b-complete semimodule. 

Note that quasisubsemimodules and A-subsemimodules may fail to be 
subsemimodules, because only the order is induced and not the correspond- 
ing addition (see Example [IS below). 

Recall that idempotent semimodules over semifields are idempotent spaces. 
In idempotent mathematics, such spaces are analogs of traditional linear 
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(vector) spaces over fields. In a similar way we use the corresponding terms 
like b-spaces, b-subspaces, b-dosed subspaces, A-subspaces, etc. 
Some examples are presented below. 

7.5. Functional semimodules. Let X be an arbitrary nonempty set and 
K be an idempotent semiring. By K{X) denote the semimodule of all 
mappings (functions) X ^ K endowed with the pointwise operations. 
By Kb{X) denote the subsemimodule of K{X) consisting of all bounded 
mappings. If i^ is a 6-complete semiring, then K{X) and K}j{X) are 6- 
complete semimodules. Note that Kb{X) is a 6-subsemimodule but not a 
6-closed subsemimodule of K{X). Given a point x G X , by 6^ denote the 
functional on K{X) that maps / to f{x). It can easily be checked that the 
functional 6x is 6- linear on K{X). 
Recall that the functional 5x is generated by the usual function 

so (p{x) = J 6x{y)(p{y)dy = sup{6x{y) Q(p{y)). Note that (5-functions form 

y 
a natural (continuous in general) basis in any typical functional semimod- 
ule. 

We say that a quasisubsemimodule of K{X) is an (idempotent) func- 
tional semimodule on the set X. An idempotent functional semimodule in 
K{X) is called b-complete if it is a 6-complete semimodule. 

A functional semimodule V C K{X) is called a functional b-semimodule 
if it is a b-subsemimodule of K{X)] a functional semimodule V C K{X) 
is called a functional A-semimodule if it is a A-subsemimodule of K{X). 

In general, a functional of the form 6^ on a functional semimodule is 
not even linear, much less 6- linear (see Example [TS] below). However, the 
following proposition holds, which is a direct consequence of our definitions. 
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Proposition 14. An arbitrary b-complete functional semimodule W on a 
set X is a b- sub semimodule of K{X) if and only if each functional of the 
form 6x (where x G X) is b-linear on W . 

Example 15. The semimodule K}){X) (consisting of all bounded mappings 
from an arbitrary set X to a 6-complete idempotent semiring K) is a func- 
tional A-semimodule. Hence it is a 6-complete semimodule over K. More- 
over, Ki){X) is a 6-subsemimodule of the semimodule K{X) consisting of 
all mappings X ^ K . 

Example 16. If X is a finite set consisting of n elements (n > 0), then 
Ki{X) = K{X) is an "n-dimensional" semimodule over K] it is denoted 
by K^. In particular, RJ^^x ^^ ^^ idempotent space over the semifield R^aa;, 
and R^ax is ^ semimodule over the semiring Rmax- Note that R^ax can 
be treated as a space over the semifield Rmax- For example, the semiring 
Rmax can be treated as a space (semimodule) over Rmax- 

Example 17. Let X be a topological space. Denote by USC{X) the set of 
all upper semicontinuous functions with values in Rmax- By definition, a 
function f{x) is upper semicontinuous if the set Xg = {x ^ X \ f{x) > s} 
is closed in X for every element s G Rmax (see, e.g., [12], Sec. 2.8). If a 
family {fa} consists of upper semicontinuous (e.g., continuous) functions 
and f{x) = inf a fa{x), then f(x) G USC{X). It is easy to check that 
USC{X) has a natural structure of an idempotent space over Rmax- More- 
over, USC{X) is a functional A-space on X and a b-space. The subspace 
USC{X)r]Kb{X) ofUSC{X) consisting of bounded (from above) functions 
has the same properties. 

Example 18. Note that an idempotent functional semimodule (and even a 
functional A-semimodule) on a set X is not necessarily a subsemimodule 
of K{X). The simplest example is the functional space (over K = Rmax) 
Conc(R) consisting of all concave functions on R with values in Rmax- 
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Recall that a function / belongs to Conc(R) if and only if the subgraph of 
this function is convex, i.e., the formula f{ax + (1 — a)y) > af{x) + (1 — 
a)f{y) is valid for < a < 1. The basic operations with G Rmax can 
be defined in an obvious way. If f^g GConc(R), then denote hy f (B g the 
sum of these functions in Conc(R). The subgraph of / © 5^ is the convex 
hull of the subgraphs of / and g. Thus f (B g does not coincide with the 
pointwise sum (i.e., max{/(a:),5'(x)}). 

Example 19. Let X be a nonempty metric space with a fixed metric r. 
Denote by Lip(X) the set of all functions defined on X with values in 
Rmax satisfying the following Lipschitz condition: 

I /(^) ifivT' 1=1 /(^) - fiy) \< r{x, y). 

where x, y are arbitrary elements of X. The set Lip(X) consists of con- 
tinuous real- valued functions (but not all of them!) and (by definition) 
the function equal to —00 = at every point x G X. The set Lip(X) 
has the structure of an idempotent space over the semifield Rmax- Spaces 
of the form Lip(X) are said to be Lipschitz spaces. These spaces are 6- 
subsemimodules in K{X). 

7.6. Integral representations of linear operators in functional semi- 
modules. Let W be an idempotent 6-complete semimodule over a 6-complete 
semiring K and V C K[X) be a 6-complete functional semimodule on X. 
A mapping A : V ^ W is called an integral operator or an operator with 
an integral representation if there exists a mapping /c : A — > VK, called the 
integral kernel (or kernel) of the operator ^4, such that 

(15) Af = sw^{f{x)Qk{x)). 

In idempotent analysis, the right-hand side of formula (11) is often written 
as J^ j{x) k{x)dx. Regarding the kernel /c, it is assumed that the set 
{f{x)Qk{x)\x G A} is bounded in W for all / G F and x G A. We denote 
the set of all functions with this property by keruywiX). In particular, if 
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W = K and ^4 is a functional, then this functional is called integral. Thus 
each integral functional can be presented in the form of a "scalar product" 
f ^ Jx /(^) ^(^) dx, where k{x) G K{X); in idempotent analysis, this 
situation is standard. 

Note that a functional of the form 6y (where y G X) is a typical integral 
functional; in this case, k{x) = 1 if x = y and k{x) = otherwise. 

We call a functional semimodule V C K{X) nondegenerate if for every 
point X G X there exists a function g E V such that g{x) = 1, and 
admissible if for every function f E V and every point x E X such that 
f{x) j^ there exists a function g eV such that g{x) = 1 and f{x)Qg :< f. 

Note that all idempotent functional semimodules over semifields are ad- 
missible (it is sufficient to set g = f{x)^^ /). 

Proposition 20. Denote by Xy the subset of X defined by the formula 
Xy = {x E X \ 3f E V : f{x) = 1}. // the semimodule V is admissible, 
then the restriction to Xy defines an embedding i : V ^ K{Xy) and its 
image i{V) is admissible and nondegenerate. 

If a mapping k : X ^ W is a kernel of a mapping A : V ^ W , then the 
mapping ky : X ^ W that is equal to k on Xy and equal to on X \ Xy 
is also a kernel of A. 

A mapping A : V ^ W is integral if and only if the mapping i-iA : 
i{A) -^ W is integral. 

In what follows, K always denotes a fixed 6-complete idempotent (basic) 
semiring. If an operator has an integral representation, this representation 
may not be unique. However, if the semimodule V is nondegenerate, then 
the set of all kernels of a fixed integral operator is bounded with respect 
to the natural order in the set of all kernels and is closed under the supre- 
mum operation applied to its arbitrary subsets. In particular, any integral 
operator defined on a nondegenerate functional semimodule has a unique 
maximal kernel. 
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An important point is that an integral operator is not necessarily b- 
linear and even linear except when V is a 6-subsemimodule of K{X) (see 
Proposition ED below) . 

If VK is a functional semimodule on a nonempty set y, then an integral 
kernel k of an operator A can be naturally identified with the function on 
X xY defined by the formula k{x^y) = {k{x)){y). This function will also 
be called an integral kernel (or kernel) of the operator A. As a result, the 
set kevny^wiX) is identified with the set kevny^wi^: Y) of all mappings k : 
X xY ^ K such that for every point x G X the mapping k^ : y \-^ k{x,y) 
lies in W and for every v G V the set {v{x) Q kx\x G X} is bounded in 
W. Accordingly, the set of all integral kernels of 6- linear operators can be 
embedded into kernv^p^(X, y). 

If V and W are functional 6-semimodules on X and y, respectively, 
then the set of all kernels of 6-linear operators can be identified with 
kerny^PT/(X, y) and the following formula holds: 

(16) A f{y) = sup {f{x)Qk{x,y))= f{x)Qk{x,y)dx. 

xex J X 

This formula coincides with the usual definition of an integral representa- 
tion of an operator. Note that formula ( jTSl ) can be rewritten in the form 

(17) Af = snp{5,{f)Qk{x)). 

Proposition 21. An arbitrary b-complete functional semimodule V on a 
nonempty set X is a functional b-semimodule on X (i.e., a b-subsemimo- 
dule of K{X)) if and only if all integral operators defined on V are b-linear. 

The following notion (definition) is especially important for our purposes. 
Let V C K{X) be a 6-complete functional semimodule over a 6-complete 
idempotent semiring K. We say that the kernel theorem holds for the semi- 
module V if every 6-linear mapping from V into an arbitrary 6-complete 
semimodule over K has an integral representation. 
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Theorem 22. Assume that a b-complete semimodule W over a b-complete 
semiring K and an admissible functional A-semimodule V C K{X) are 
given. Then every b-linear operator A : V ^ W has an integral represen- 
tation of the form (TTHI). In particular, ifW is afunctional b-semimodule on 
a setY , then the operator A has an integral representation of the form (TTBl). 
Thus for the semimodule V the kernel theorem holds. 

Remark 23. Examples of admissible functional A-semimodules (and A- 
spaces) appearing in Theorem [22] are presented above, see, e.g., exam- 
ples [15] - [El Thus for these functional semimodules and spaces V over K, 
the kernel theorem holds and every 6-linear mapping V into an arbitrary 
6-complete semimodule W over K has an integral representation (1161). Re- 
call that every functional space over a 6-complete semifield is admissible, 
see above. 

7.7. Nuclear operators and their integral representations. Let us 

introduce some important definitions. Assume that V and W are 6-complete 
semimodules. A mapping g : V ^ W is called one- dimensional (or a map- 
ping of rank 1) if it is of the form v i— > 0('u) w;, where is a 6-linear 
functional on V and w G W. A mapping g is called b-nuclear if it is 
the sum (i.e., supremum) of a bounded set of one-dimensional mappings. 
Since every one-dimensional mapping is 6-linear (because the functional 
is 6-linear), every b-nuclear operator is b-linear (see Corollary [T2] above). 
Of course, 6-nuclear mappings are closely related to tensor products of 
idempotent semimodules, see pS] . 

By 00w; we denote the one-dimensional operator v i— > (j){v) Qw. In fact, 
this is an element of the corresponding tensor product. 

Proposition 24. The composition (product) of a b-nuclear and a b-linear 
mapping or of a b-linear and a b-nuclear mapping is a b-nuclear operator. 

Theorem 25. Assume that W is a b-complete semimodule over a b-complete 
semiring K and V C K{X) is a functional b-semimodule. If every b-linear 
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functional on V is integral, then a b-linear operator A : V ^ W has an 
integral representation if and only if it is b-nuclear. 

7.8. The 6- approximation property and 6-nuclear semimodules and 
spaces. We say that a 6-complete semimodule V has the b- approximation 
property if the identity operator id:y — > V is 6-nuclear (for a treatment 
of the approximation property for locally convex spaces in the traditional 
functional analysis, see [86j). 

Let V be an arbitrary 6-complete semimodule over a 6-complete idempo- 
tent semiring K. We call this semimodule a b-nuclear semimodule if any 
6-linear mapping of V to an arbitrary 6-complete semimodule W over K is 
a 6-nuclear operator. Recall that, in the traditional functional analysis, a 
locally convex space is nuclear if and only if all continuous linear mappings 
of this space to any Banach space are nuclear operators, see [HB] . 



Proposition 26. Let V be an arbitrary b-complete semimodule over a b- 
complete semiring K . The following statements are equivalent: 

1 the semimodule V has the b- approximation property; 

2 every b-linear mapping from V to an arbitrary b-complete semimodule 
W over K is b-nuclear; 

3 every b-linear mapping from an arbitrary b-complete semimodule W 
over K to the semimodule V is b-nuclear. 

Corollary 27. An arbitrary b-complete semimodule over a b-complete semir- 
ing K is b-nuclear if and only if this semimodule has the b- approximation 
property. 

Recall that, in the traditional functional analysis, any nuclear space has 
the approximation property but the converse is not true. 

Concrete examples of 6-nuclear spaces and semimodules are described 
in Examples [151 [IS] and [IS] (see above). Important 6-nuclear spaces and 
semimodules (e.g., the so-called Lipschitz spaces and semi-Lipschitz semi- 
modules) are described in [57]. In this paper there is a description of 
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all functional 6-semimodules for which the kernel theorem holds (as semi- 
Lipschitz semimodules); this result is due to G. B. Shpiz. 

It is easy to show that the idempotent spaces USC{X) and Conc(R 
(see Examples [IT] and [IE]) are not 6-nuclear (however, for these spaces the 
kernel theorem is true). The reason is that these spaces are not functional 
6-spaces and the corresponding (5-functionals are not 6-linear (and even 
linear). 

7.9. Kernel theorems for functional 6-seniiniodules. Let V C K{X) 

be a 6-complete functional semimodule over a 6-complete semiring K . Re- 
call that for V the kernel theorem holds if every 6-linear mapping of this 
semimodule to an arbitrary 6-complete semimodule over K has an integral 
representation. 

Theorem 28. Assume that a b-complete semiring K and a nonempty set 
X are given. The kernel theorem holds for any functional b-semimodule 

V C K{X) if and only if every b-linear functional on V is integral and the 
semimodule V is b-nuclear, i.e., has the b- approximation property. 

Corollary 29. If for a functional b-semimodule the kernel theorem holds, 
then this semimodule is b-nuclear. 

Note that the possibility to obtain an integral representation of a func- 
tional means that one can decompose it into a sum of functionals of the 
form 5x. 

Corollary 30. Assume that a b-complete semiring K and a nonempty 
set X are given. The kernel theorem holds for a functional b-semimodule 

V C K{X) if and only if the identity operator id: V ^V is integral. 

7.10. Integral representations of operators in abstract idempo- 
tent semimodules. In this subsection, we examine the following problem: 
when a 6-complete idempotent semimodule V over a 6-complete semiring is 
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isomorphic to a functional 5-semimodule W such that the kernel theorem 
holds for W. 

Assume that F is a 6-complete idempotent semimodule over a 6-complete 
semiring K and is a 6-linear functional defined on V. We call this func- 
tional a 6-functional if there exists an element v eV such that 

(j){w) Q V ^ w 

for every element w E V. It is easy to see that every functional of the form 
6x is a ^-functional in this sense (but the converse is not true in general). 
Denote by A(V) the set of all ^- functional on V. Denote by ^a the 
natural mapping V — > K{A{V)) defined by the formula 

for all (f) G A(y). We say that an element -u G F is pointlike if there exists 
a 6-linear functional such that 4>{w) Q v ^ w for all w E V. The set of 
all pointlike elements of V will be denoted by P{V). Recall that hj (f)Q v 
we denote the one-dimensional operator w i— > (f){w) v. 

The following assertion is an obvious consequence of our definitions (in- 
cluding the definition of the standard order) and the idempotency of our 
addition. 

Remark 31. If a one-dimensional operator (f)Qv appears in the decomposi- 
tion of the identity operator on V into a sum of one-dimensional operators, 
then (f) E A{V) and v E P{V). 

Denote by id and Id the identity operators on V and ^a(^), respectively. 

Proposition 32. // the operator id is b-nuclear, then i^ is an embedding 
and the operator Id is integral. 

If the operator i/^ is an embedding and the operator Id is integral, then 
the operator id is b-nuclear. 

Theorem 33. A b-complete idempotent semimodule V over a b-complete 
idempotent semiring K is isomorphic to a functional b-semimodule for 
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which the kernel theorem holds if and only if the identity mapping on V is 
a b-nuclear operator, i.e., V is a b-nuclear semimodule. 

The following proposition shows that, in a certain sense, the embedding 
za is a universal representation of a 6-nuclear semimodule in the form of a 
functional 6-semimodule for which the kernel theorem holds. 

Proposition 34. Let K be a b-complete idempotent semiring, X be a 
nonempty set, and V C K{X) be a functional b-semimodule on X for 
which the kernel theorem holds. Then there exists a natural mapping 
i : X ^ /S.{V) such that the corresponding mapping i^ : K{A{V)) — > K{X) 
is an isomorphism of iA{V) onto V. 

8. The dequantization transform, convex geometry and the 

Newton polytopes 

Let X be a topological space. For functions f{x) defined on X we shall 
say that a certain property is valid almost everywhere (a.e.) if it is valid 
for all elements x of an open dense subset of X. Suppose X is C^ or R*^; 
denote by R" the set x = { (xi, . . . , Xn) ^ X \ Xi > for i = 1,2, . . . ,n. 
For X = {xi, . . . ,Xn) G X we set exp(a:) = (exp(a:i), . . . ,exp(a:„)); so if 
X ell^, then exp(a:) G R" . 

Denote by T{C^) the set of all functions defined and continuous on an 
open dense subset U C C"' such that U D R*^. It is clear that T{C^) 
is a ring (and an algebra over C) with respect to the usual addition and 
multiplications of functions. 

For / G J-'(C'^) let us define the function fh by the following formula: 

(18) fh{x) = hlog\f{exp{x/h))\, 

where /i is a (small) real positive parameter and x G R". Set 

(19) />)= hm A(x), 

if the right-hand side of (1191 ) exists almost everywhere. 
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We shall say that the function f{x) is a dequantization of the function 
f{x) and the map f{x) i— > f{x) is a dequantization transform. By con- 
struction, fh{x) and f{x) can be treated as functions taking their values 
in Rmax- Note that in fact fh{x) and f{x) depend on the restriction of 
/ to R" only; so in fact the dequantization transform is constructed for 
functions defined on R" only. It is clear that the dequantization transform 
is generated by the Maslov dequantization and the map a: i— > \x\. 

Of course, similar definitions can be given for functions defined on R"^ 
and IU\_. If s = 1/h, then we have the following version of ( TTHI ) and (HM: 

(20) f{x)= lim(l/s)log|/(e^^)|. 

Denote by df the subdifferential of the function / at the origin. 
If / is a polynomial we have 

df = {ve R" I {v, x) < f{x) Vx G R" }. 

It is well known that all the convex compact subsets in R"^ form an idem- 
potent semiring S with respect to the Minkowski operations: for a, /5 G <S 
the sum a /3 is the convex hull of the union a U /3; the product a /3 is 
defined in the following way: a0/5 = {a:|a: = a + 6, where a G a^b G jS, 
see Fig El In fact S is an idempotent linear space over Rmax- 

Of course, the Newton polytopes of polynomials in n variables form a 
subsemiring J\f in S. If /, g are polynomials, then d{fg) = df Q d^; 
moreover, if / and g are "in general position", then d{f + g) = df ® 
dg. For the semiring of all polynomials with nonnegative coefficients the 
dequantization transform is a homomorphism of this "traditional" semiring 
to the idempotent semiring TV. 
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Theorem 35. If f is a polynomial, then the subdijjerential d f of f at the 
origin coincides with the Newton polytope of f. For the semiring of poly- 
nomials with nonnegative coefficients, the transform f \-^ df is a homo- 
morphism of this semiring to the semiring of convex polytopes with respect 
to the Minkowski operations (see above). 

Using the dequantization transform it is possible to generahze this resuh 
to a wide class of functions and convex sets, see below and [55]. 

8.1. Dequantization transform: algebraic properties. Denote by V 
the set R"^ treated as a linear Euclidean space (with the scalar product 

(x, y) = xiyi + X2y2 + • • • + Xny-n) and set y+ = R" . We shall say that 
a function / G J^{C^) is dequantizahle whenever its dequantization f{x) 
exists (and is defined on an open dense subset of V). By X>(C") denote 
the set of all dequantizahle functions and by ViV) denote the set { / | 
/ G P(C") }. Recall that functions from P(C") (and V{V)) are defined 
almost everywhere and f = g means that f{x) = g{x) a.e., i.e., for x 
ranging over an open dense subset of C" (resp., of V). Denote by 'D+(C'^) 
the set of all functions / G V{C"') such that f{xi, . . . ,x„) > if x^ > 
for ? = 1, . . . , n; so / G X>+(C") if the restriction of f to V+ = R*^ is a 
nonnegative function. By P+(y) denote the image of X>+(C"^) under the 
dequantization transform. We shall say that functions /, (/ G P(C") are in 
general position whenever f{x) ^ g{x) for x running an open dense subset 

of y. 




Figure 3. Algebra of convex subsets. 
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Theorem 36. For functions f^g ^ I>(C") and any nonzero constant c, 
the following equations are valid: 

2) |/H/;cf = /;^=0; 

3) {f + g){x) = max{/(a:),^(x)} a.e. if f and g are nonnegative on V+ 
(i.e., f^gE Vj^{C'^)) or f and g are in general position. 

Left-hand sides of these equations are well-defined automatically. 



Corollary 37. The set P+(C'^) has a natural structure of a semiring with 
respect to the usual addition and multiplication of functions taking their 
values in C. The set V+{V) has a natural structure of an idempotent 
semiring with respect to the operations (/ g){x) = m.ax{f{x)^g{x)}, 
{f & g){x) = f{x) + g{x); elements of T>_^_{V) can be naturally treated 
as functions taking their values in Rmax- The dequantization transform 
generates a homomorphism from T>j^{C^) to T>+{V). 

8.2. Generalized polynomials and simple functions. For any nonzero 

number a E C and any vector d = (di, . . . , dn) eV = R"^ we set rUa^dix) = 
aYYi^i X-'; functions of this kind we shall call generalized monomials. Gen- 
eralized monomials are defined a.e. on C" and on y+, but not on V unless 
the numbers di take integer or suitable rational values. We shall say that a 
function / is a generalized polynomial whenever it is a finite sum of linearly 
independent generalized monomials. For instance, Laurent polynomials 
and Puiseax polynomials are examples of generalized polynomials. 

As usual, for x^y E V we set (x, y) = xiyi + • • • + Xny-n- The following 
proposition is a result of a trivial calculation. 

Proposition 38. For any nonzero number a E V = C and any vector 
d eV = R"^ we have {fn^d)h{x) = (<i, x) + h log \a\ . 

Corollary 39. /// is a generalized monomial, then f is a linear function. 
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Recall that a real function j9 defined on F = R" is sublinearifp = sup^ j9q,, 
where {pa} is a collection of linear functions. Sublinear functions defined 
everywhere on y = R"^ are convex; thus these functions are continuous, 
see jH] • We discuss sublinear functions of this kind only. Suppose ^ is a 
continuous function defined on V, then p is sublinear whenever 

1) p{x + y) < p{x) + p{y) for all x^y ^V] 

2) p{cx) = cp{x) for all x G y, c G R+. 

So if j3i, P2 are sublinear functions, then pi + j92 is a sublinear function. 

We shall say that a function / G J^(C") is simple^ if its dequantization / 
exists and a.e. coincides with a sublinear function; by misuse of language, 
we shall denote this (uniquely defined everywhere on V) sublinear function 
by the same symbol /. 

Recall that simple functions / and g are in general position if f{x) ^ 'g{x) 
for all X belonging to an open dense subset of V. In particular, generalized 
monomials are in general position whenever they are linearly independent. 

Denote by SimiCT) the set of all simple functions defined on V and de- 
note by 5'zm+(C") the set Sim{C'') nV+{C''). By Shl{V) denote the set 
of all (continuous) subhnear functions defined on "1/ = R" and by Shlj^{V) 
denote the image 5'zm+(C") of Sim+{C'^) under the dequantization trans- 
form. 

The following statements can be easily deduced from Theorem 8.2 and 
definitions. 

Corollary 40. The set Sim j^{C") is a suhsemiring ofV+(C^) and Sbl+{V) 
is an idempotent suhsemiring o/P+ ("!/). The dequantization transform 
generates an epimorphism of Sim^{C^) onto Sbl+{V). The set Sbl{V) 
is an idempotent semiring with respect to the operations (/ © g){x) = 
miix{f{x),g{x)}, {fQg){x) = f{x) + g{x). 

Corollary 41. Polynomials and generalized polynomials are simple func- 
tions. 
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We shall say that functions f^g ^ ^(^) ^^^ asymptotically equivalent 
whenever f = g\ any simple function / is an asymptotic monomial when- 
ever / is a linear function. A simple function / will be called an asymptotic 
polynomial whenever / is a sum of a finite collection of nonequivalent as- 
ymptotic monomials. 

Corollary 42. Every asymptotic polynomial is a simple function. 

Example 43. Generalized polynomials, logarithmic functions of (general- 
ized) polynomials, and products of polynomials and logarithmic functions 
are asymptotic polynomials. This follows from our definitions and for- 
mula (US). 

8.3. SubdifFerentials of sublinear functions. We shall use some ele- 
mentary results from convex analysis. These results can be found, e.g., 
in |M], ch. 1, §1. 
For any function p E Sbl{V) we set 

(21) dp={v eV \{v,x) <p{x)yx eV}. 

It is well known from convex analysis that for any sublinear function p 
the set dp is exactly the subdifferential of p at the origin. The following 
propositions are also known in convex analysis. 

Proposition 44. Suppose pi,p2 G Sbl{V), then 

1) d{pi+p2) = dpiQdp2 = {v gV \ V = 'Ui+'U2, where vi G dpi^ V2 G dp2 }; 

2) d{m.ax{pi{x) , p2{x)}) = dpi © dp2. 

Recall that dpi © dp2 is a convex hull of the set dpi U dp2. 

Proposition 45. Suppose p G Sbl{V). Then dp is a nonempty convex 
compact subset ofV. 

Corollary 46. The map p \-^ dp is a homomorphism of the idempotent 
semiring Sbl{V) (see Corollary f^^ to the idempotent semiring S of all 



convex compact subsets ofV (see Subsection \8A\ above). 
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8.4. Newton sets for simple functions. For any simple function / G 
SimiC) let us denote by N{f) the set d{f). We shall caU N{f) the Newton 
set of the function /. 

Proposition 47. For any simple function f, its Newton set N{f) is a 
nonempty convex compact subset ofV. 

This proposition follows from Proposition US and definitions. 

Theorem 48. Suppose that f and g are simple functions. Then 

1) N{fg) = N{f)QN{g) = {veV\v = V1+V2 with vi G N{f), V2 G N{g) }; 

2) A^(/ + g) = N{f) N{g), if /i and /2 are in general position or 
/i,/2 G Sim^{C^) {recall that N{f) © N{g) is the convex hull of 
Nif)UN{g)). 

This theorem follows from Theorem [3Sl Proposition SH and definitions. 

Corollary 49. The map f 1— > N{f) generates a homomorphism from 
Sim+iC") to S. 

Proposition 50. Let f = ma^d{x) = oHiLi^i' ^^ ^ monomial; here d = 
(di, . . . , dn) G y = R" and a is a nonzero complex number. Then N{f) = 

{d}. 

This follows from Proposition EHl Corollary [33 and definitions. 

Corollary 51. Let f = YlideD'^^aa-.d be a polynomial. Then N{f) is the 
polytope (Bdeoid}, i.e. the convex hull of the finite set D. 

This statement follows from Theorem SHI and Proposition [5U1 Thus in this 
case N{f) is the well-known classical Newton polytope of the polynomial 

Now the following corollary is obvious. 

Corollary 52. Let f be a generalized or asymptotic polynomial. Then its 
Newton set N{f) is a convex polytope. 
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Example 53. . Consider the one dimensional case, i.e., V = 11 and suppose 

/i = a^x" + an-ix''-'^ H h ao and /2 = 6^3:™ + 6m-i3:"'"^ H h 6o, 

where a„ 7^ 0, 6^ 7^ 0, ao 7^ 0, 60 7^ 0. Then N{fi) is the segment [0, n] 
and N{f2) is the segment [0,m]. So the map / 1— > A^(/) corresponds to 
the map / 1— > deg(/), where deg(/) is a degree of the polynomial /. In 
this case Theorem 2 means that deg{fg) = deg / + deg g and deg(/ + (/) = 
maxjdeg /, deg^^} = max{n, m} if a^ > 0, 6^ > or / and g are in general 
position. 

9. Dequantization of set functions and measures on metric 

SPACES 
The following results are presented in [5^ . 



Example 54. Let M be a metric space, S its arbitrary subset with a compact 
closure. It is '' 
p has volume 



closure. It is well-known that a Euclidean d-dimensional ball Bp of radius 



r(i/2)^ 
'^^ ~ r(i + d/2)' 



yoid{Bp) = ^,; ; ^^, p' 



d 



where (i is a natural parameter. By means of this formula it is possible to 
define a volume of Bp for any real d. Cover S* by a finite number of balls 
of radii pm. Set 

Vd{S) := lim inf ^Yo\d{BpJ. 

m 

Then there exists a number D such that Vd{S) = for d > D and ^^(5') = 
00 for d < D. This number D is called the Hausdorff-Besicovich dimension 
(or HB-dimension) of 5, see, e.g., [67]. Note that a set of non- integral HB- 
dimension is called a fractal in the sense of B. Mandelbrot. 

Theorem 55. Denote by Afp{S) the minimal number of balls of radius p 
covering S . Then 

D{S) = hm \ogp{Mp{Sr'), 
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where D{S) is the HB-dimension of S. Set p = e^^ , then 

D{S)= iim {l/s)-logKM-s)iS). 

So the HB-dimension D{S) can be treated as a result of a dequantization 
of the set function Mp{S). 

Example 56. Let /i be a set function on M (e.g., a probability measure) 
and suppose that p{Bp) < oo for every ball Bp. Let B^^p be a ball of radius 
p having the point a: G M as its center. Then define Px{p) '•= /^(5x,p) and 
let p = e^^ and 

D,^,:= hm -(l/s)-log(|/i,(e--^)|). 



-CXD 



This number could be treated as a dimension of M at the point x with 
respect to the set function p. So this dimension is a result of a dequanti- 
zation of the function Px{p)-i where x is fixed. There are many dequantiza- 
tion procedures of this type in different mathematical areas. In particular, 
V.P. Maslov's negative dimension (see [S7]) can be treated similarly. 

10. Dequantization of geometry 

An idempotent version of real algebraic geometry was discovered in the 
report of O. Viro for the Barcelona Congress jSl]- Starting from the idem- 
potent correspondence principle O. Viro constructed a piecewise-linear ge- 
ometry of polyhedra of a special kind in finite dimensional Euclidean spaces 
as a result of the Maslov dequantization of real algebraic geometry. He indi- 
cated important applications in real algebraic geometry (e.g., in the frame- 
work of Hubert's 16th problem for constructing real algebraic varieties with 
prescribed properties and parameters) and relations to complex algebraic 
geometry and amoebas in the sense of L M. Gelfand, M. M. Kapranov, 
and A. V. Zelevinsky, see [2SIE5]- Then complex algebraic geometry was 
dequantized by G. Mikhalkin and the result turned out to be the same; this 
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new 'idempotent' (or asymptotic) geometry is now often called the tropical 
algebraic geometry, see, e.g., P2lli5lli^[53l[721[75] . 

There is a natural relation between the Maslov dequantization and amoe- 
bas. 

Suppose (C*)" is a complex torus, where C* = C\{0} is the group of 
nonzero complex numbers under multiplication. For z = (21,..., z„) G 
(C*)'^ and a positive real number h denote by Logf^{z) = /ilog(|2;|) the 
element 

{hlog \zi\,hlog \z2\: . . . , /ilog l^^l) e R". 

Suppose V C (C*)" is a complex algebraic variety; denote by Ah{V) the set 
Logj^{V). If h = 1, then the set A{V) = Ai{V) is called the amoeba of V; 
the amoeba A{V) is a closed subset of R" with a non-empty complement. 
Note that this construction depends on our coordinate system. 

For the sake of simplicity suppose F is a hypersurface in (C*)'^ defined 
by a polynomial /; then there is a deformation h \-^ fh of. this polynomial 
generated by the Maslov dequantization and fh = f for h = 1. Let "14 C 
(C*)'* be the zero set of fh and set Ah{Vh) = Log/j("\4). Then there exists a 
tropical variety Tro{V) such that the subsets Ah{Vh) C R" tend to Tro{V) 
in the Hausdorff metric as /i — > 0. The tropical variety Tro{V) is a result 
of a deformation of the amoeba ^(^) and the Maslov dequantization of 
the variety V. The set Tro{y) is called the skeleton of A{V). 





(a) (b) (c) 

Figure 4. Tropical line and deformations of an amoeba 
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Example 57. For the line V = {{x,y) ^ (C*)^ | x+y+1 = } the piecewise- 
hnear graph Tro{V) is a tropical line, see Fig. [11(a). The amoeba A{V) is 
represented in Fig. [11(b), while Fig. 111(c) demonstrates the corresponding 
deformation of the amoeba. 

11. Some semiring constructions and the matrix Bellman 

equation 

11.1. Complete idempotent semirings and examples. Recall that a 
partially ordered set S is complete if for every subset T C S there exist 
elements sup T G S and inf T G S. We say that an idempotent semiring S 
is complete if it is complete as an ordered set with respect to the standard 
order. Of course, any a-complete semiring (see subsect. 17.11 ) is complete. 
The most well-known and important examples are "numerical semirings" 
consisting of (a subset of) real numbers and ordered by the usual linear 
order <. 

Example 58. Consider the semiring Rmax = R-max U {oo} with standard 
operations © = max, = + and neutral elements = — oo, 1 = 0, 
X < oo, X © oo = oo for all X, X oo = oo a: = oo if a: 7^ 0, and 

oo = oo 0. The semiring Rmax is complete and a-complete. The 
semiring Rmin = Rmin U {— oo} with obvious operations is also complete; 
Rmin and Rmax are isomorphic. 

Example 59. Consider the semiring 5'||Jaxmin defined on the real interval 
[a, b] with operations © = max, = min and neutral elements = a and 

1 = 6. The semiring is complete and a-complete. Set ^max,!!!!!! = 'S'^^^ax min 
with a = — oo and b = +oo. If — oo < a < b < +oo then »S',^'i„,;„ and 
5'max,min are Isomorphlc. 

Example 60. The Boolean algebra 5 = {0,l}isa complete and a-complete 
semifield consisting of two elements. 
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11.2. Closure operations. Let a semiring S be endowed with a partial 
unary closure (or Kleene) operation * such that x ^ y implies x* :< y* and 
X* = 1 © {x* x) = 1 © (x X*) on its domain of definition. In particular, 
0* = 1 by definition. These axioms imply that x* = l©x©x^©---© 
{x* x^) if n ^ 1. Thus x* can be considered as a 'regularized sum' of 
the series x* = 1 © x © x^ ® . . . ; in an idempotent semiring, by definition, 
X* = sup{l,a:,x^, . . . } if this supremum exists. So if S is complete, then 
the closure operation is well-defined for every element x E S. 

In numerical semirings the operation * is defined as follows: x* = (1 — 
x)^^ if x -< 1 in R+, or R+ and x* = oo if x ^ 1 in R+; x* = 1 if x ^ 1 
in Rniax and Rmax, X* = oo if X ^ 1 in R^ax, x* = 1 for all x in S!^'^l^-^. 
In all other cases x* is undefined. Note that the closure operation is very 
easy to implement. 

11.3. Matrices over semirings. Denote by MaimniS) a set of all matri- 
ces A = [ttij) with m rows and n columns whose coefficients belong to a 
semiring S. The sum ^4©^ of matrices A, B E MdutmniS) and the product 
AB of matrices A G MatimiS) and B G Mdutmn{S) are defined according 
to the usual rules of linear algebra: A^ B = (aij © bij) G MaimniS) and 

where A G MaiimiS) and B G Mai^nniS). Note that we write AB instead 
ofAQB. 

If the semiring S is ordered, then the set Mdutmn{S) is ordered by the 
relation A = [aij) :< B = [bij) iff aij :< bij in S for all 1 ^ z ^ m, 

1 ^ j ^ ^• 

The matrix multiplication is consistent with the order ^ in the following 
sense: ff A, A' G Mat;^(5'), B,B' G Mat^„(5') and A ^ A\ B ^ B\ then 
AB :< A'B' in Mdutin{S). The set Matn„(S') of square (n x n) matrices over 
an idempotent semiring S forms a idempotent semiring with a zero element 
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O = {oij)^ where Oij = 0, 1 ^ ^,j ^ n, and a unit element / = (%), where 
5ij = I'lii = i and 5ij = otherwise. 

The set Mat^„ is an example of a noncommutative semiring if n > 1. 

The closure operation in matrix semirings over an idempotent semiring S 
can be defined inductively (another way to do that see in [22] and below): 
A* = (ail)* = (fi^ii) in Matii(5') and for any integer n > 1 and any matrix 

^ ^ (All Au 

\A21 A22^ 

where An e Matfcfc(5'), Au G Matkn-k{S), A21 G Md.tn-kk{S), A22 G 
Matn-kn-kiS), 1 ^ /c ^ n, by defintion, 

/Al, e Al,Ai2D*A2iAl, Al.AuD' 

(22) A* = 

\ D*A2iA*n D* 

where D = A22 © ^21^11^12- It can be proved that this definition of 
A* implies that the equality A* = A* A © / is satisfied and thus A* is a 
'regularized sum' of the series / © ^4 © ^4^ © . . . . 

Note that this recurrence relation coincides with the formulas of escalator 
method of matrix inversion in the traditional linear algebra over the field of 
real or complex numbers, up to the algebraic operations used. Hence this 
algorithm of matrix closure requires a polynomial number of operations 
in n. 

11.4. Discrete stationary Bellman equations. Let 5 be a semiring. 
The discrete stationary Bellman equation has the form 

(23) X = AXeB, 

where A G Mat„„(5'), X, B E Mat„5(5'), and the matrix X is unknown. Let 
A* be the closure of the matrix A. It follows from the identity A* = 74*^4©/ 
that the matrix A^'B satisfies this equation; moreover, it can be proved that 
for idempotent semirings this solution is the least in the set of solutions to 
equation ([23]) with respect to the partial order in Mat„5(5'). 
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Equation ( 1231 ) over max-plus semiring arises in connection with Bellman 
optimality principle and discretization of Hamilton-Jacobi equations, see 
e.g., [ZI]. It is also intimately related with optimization problems on graphs 



to be discussed below. 

11.5. Weighted directed graphs and matrices over semirings. Sup- 
pose that 5* is a semiring with zero and unity 1. It is well-known that any 
square matrix A = [aij) G Mat„„(5') specifies a weighted directed graph. 
This geometrical construction includes three kinds of objects: the set X 
of n elements xi, . . . , a:„ called nodes, the set F of all ordered pairs (xj, Xj) 
such that ttij 7^ called arcs, and the mapping A: T ^ S such that 
A{xi,Xj) = aij. The elements aij of the semiring S are called weights of 
the arcs. See Fig. [5] 

Conversely, any given weighted directed graph with n nodes specifies a 
unique matrix A G Mat„r7,(»S'). 

This definition allows for some pairs of nodes to be disconnected if the 
corresponding element of the matrix ^4 is and for some channels to be 
"loops" with coincident ends if the matrix A has nonzero diagonal elements. 




Figure 5. A weighted directed graph. 
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This concept is convenient for analysis of parallel and distributed computa- 
tions and design of computing media and networks (see, e.g., [51I151ESIE7] ) . 
Recall that a sequence of nodes of the form 

P= (?/o,2/i,---,2/fc) 

with k ^ and (?/i, yi+i) G F, ? = 0, . . . , /c — 1, is called a path of length k 
connecting ijq with i/k. Denote the set of all such paths by Pk{yo, Vk)- The 
weight A{p) of a path p G Pk{yo: Vk) is defined to be the product of weights 
of arcs connecting consecutive nodes of the path: 

A{p) = A{yo, ?/i) • • • A{yk^i, yk). 

By definition, for a 'path' p G PQ{xi^Xj) of length A; = the weight is 1 if 
i = i and otherwise. 

For each matrix A G Mat,^^(5') define ^4^ = / = {5ij) (where 5ij = 1 if 
i = j and 5ij = otherwise) and A^ = AA^^^, k ^ I. Let a-- be the 
(i,j)th element of the matrix A^. It is easily checked that 

^ij ~ \±y Ojioh ■ • ■ Ciik-iik- 

io^i,ik=j 
l!^ii,...,ik-i^n 



(k) 

hj 
of length k connecting the node Xi„ = Xi with Xi^ = Xj. 



Thus a- J is the supremum of the set of weights corresponding to all paths 



Denote the elements of the matrix A* by a--\, z, j = 1, . . . , n; then 



«^u = 



-'J- 

I ■ ■ _ 

© © A{p). 

The closure matrix A* solves the well-known algebraic path problem, 
which is formulated as follows: for each pair [xi^Xj) calculate the supre- 
mum of weights of all paths (of arbitrary length) connecting node Xi with 
node Xj. The closure operation in matrix semirings has been studied ex- 
tensively (see, e.g., |Il[2l[Ml[Ill[I5l[2DH22l[2SH3DlE3l[31,[5Sl and references 
therein). 
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Example 61 (The shortest path problem.). Let S = Rmin, so the weights 
are real numbers. In this case 

A{p) = A{yo, yi) + A{yi, 7/2) H h A{yk-i, yk)- 

If the element Ojj specifies the length of the arc (xj, Xj) in some metric, 
then a\- is the length of the shortest path connecting Xi with Xj. 

Example 62 (The maximal path width problem.). Let ^ = RU {0, 1} with 
= max, = min. Then 

4?^= , pax A{p), A{p) = min{A{yQ,yi),...,A{yk^i,yk)). 
pe U Pk(xi,Xj) 

If the element aij specifies the "width" of the arc (xi,Xj), then the width 
of a path p is defined as the minimal width of its constituting arcs and the 
element a-- gives the supremum of possible widths of all paths connecting 
Xi with Xj. 

Example 63 (A simple dynamic programming problem.). Let S = Rmax 
and suppose aij gives the profit corresponding to the transition from Xi to 
Xj. Define the vector B = (bi) G Mat„i(Rmax) whose element bi gives the 
terminal profit corresponding to exiting from the graph through the node 
Xi. Of course, negative profits (or, rather, losses) are allowed. Let m be 
the total profit corresponding to a path p G Pk{xi,Xj), i.e. 

m = A{p) + bj. 

Then it is easy to check that the supremum of profits that can be achieved 
on paths of length k beginning at the node Xi is equal to {A^B)i and the 
supremum of profits achievable without a restriction on the length of a 
path equals {A*B)i. 

Example 64 (The matrix inversion problem.). Note that in the formulas 
of this section we are using distributivity of the multiplication with 
respect to the addition but do not use the idempotency axiom. Thus 
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the algebraic path problem can be posed for a nonidempotent semiring S 
as well (see, e.g., [Bl])- For instance, if 5 = R, then 

A* = I + A + A'^ + --- = {I- A)-\ 

If 1 1^4 1 1 > 1 but the matrix / — ^4 is invertible, then this expression defines 
a regularized sum of the divergent matrix power series X^j>o^*- 

There are many other important examples of problems (in different areas) 
related to algorithms of linear algebra over semirings (transitive closures 
of relations, accessible sets, critical paths, paths of greatest capacities, the 
most rehable paths, interval and other problems), see pi[21[5ti71 [T^[nHT7i 



We emphasize that this connection between the matrix closure operation 
and solution to the Bellman equation gives rise to a number of different 
algorithms for numerical calculation of the closure matrix. All these al- 
gorithms are adaptations of the well-known algorithms of the traditional 
computational linear algebra, such as the Gauss-Jordan elimination, var- 
ious iterative and escalator schemes, etc. This is a special case of the 
idempotent superposition principle. 

In fact, the theory of the discrete stationary Bellman equation can be 
developed using the identity A* = AA* © / as an additional axiom without 
any substantial interpretation (the so-called closed semirings, see, e.g., [71, 

mmm)- 

12. Universal algorithms 

Computational algorithms are constructed on the basis of certain prim- 
itive operations. These operations manipulate data that describe "num- 
bers." These "numbers" are elements of a "numerical domain," i.e., a 
mathematical object such as the field of real numbers, the ring of integers, 
or an idempotent semiring of numbers. 
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In practice elements of the numerical domains are replaced by their com- 
puter representations, i.e., by elements of certain finite models of these 
domains. Examples of models that can be conveniently used for com- 
puter representation of real numbers are provided by various modifications 
of floating point arithmetics, approximate arithmetics of rational num- 
bers [52], and interval arithmetics. The difference between mathematical 
objects ("ideal" numbers) and their finite models (computer representa- 
tions) results in computational (e.g., rounding) errors. 

An algorithm is called universal if it is independent of a particular 
numerical domain and/or its computer representation. A typical exam- 
ple of a universal algorithm is the computation of the scalar product 
(x, y) of two vectors x = (xi, . . . , Xn) and y = (?/i, . . . , yn) by the for- 
mula (x, y) = xiyi + • • • + Xnyn- This algorithm (formula) is independent 
of a particular domain and its computer implementation, since the formula 
is well-defined for any semiring. It is clear that one algorithm can be more 
universal than another. For example, the simplest Newton-Cotes formula, 
the rectangular rule, provides the most universal algorithm for numerical 
integration; indeed, this formula is valid even for idempotent integration 
(over any idempotent semiring, see above and |5l [55ll5UlH0lll2] - liil[5Tl[52] - IB5] . 



Other quadrature formulas (e.g., combined trapezoid rule or the Simpson 
formula) are independent of computer arithmetics and can be used (e.g., 
in an iterative form) for computations with arbitrary accuracy. In con- 
trast, algorithms based on Gauss- Jacobi formulas are designed for fixed 
accuracy computations: they include constants (coefficients and nodes of 
these formulas) defined with fixed accuracy. Certainly, algorithms of this 
type can be made more universal by including procedures for computing 
the constants; however, this results in an unjustified complication of the 
algorithms. 

Computer algebra algorithms used in such systems as Mathematica, 
Maple, REDUCE, and others are highly universal. Most of the standard 
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algorithms used in linear algebra can be rewritten in such a way that they 
will be valid over any field and complete idempotent semiring (including 
semirings of intervals; see below and [ J5S1[5UII^ , where an interval version of 
the idempotent linear algebra and the corresponding universal algorithms 
are discussed). 

As a rule, iterative algorithms (beginning with the successive approxi- 
mation method) for solving differential equations (e.g., methods of Euler, 
Euler-Cauchy, Runge-Kutta, Adams, a number of important versions of 
the difference approximation method, and the like), methods for calcu- 
lating elementary and some special functions based on the expansion in 
Taylor's series and continuous fractions (Pade approximations) and others 
are independent of the computer representation of numbers. 

Calculations on computers usually are based on a floating-point arith- 
metic with a mantissa of a fixed length; i.e., computations are performed 
with fixed accuracy. Broadly speaking, with this approach only the rela- 
tive rounding error is fixed, which can lead to a drastic loss of accuracy 
and invalid results (e.g., when summing series and subtracting close num- 
bers). On the other hand, this approach provides rather high speed of 
computations. Many important numerical algorithms are designed to use 
floating-point arithmetic (with fixed accuracy) and ensure the maximum 
computation speed. However, these algorithms are not universal. The 
above mentioned Gauss-Jacobi quadrature formulas, computation of ele- 
mentary and special functions on the basis of the best polynomial or ratio- 
nal approximations or Pade-Chebyshev approximations, and some others 
belong to this type. Such algorithms use nontrivial constants specified with 
fixed accuracy. 

Recently, problems of accuracy, reliability, and authenticity of computa- 
tions (including the effect of rounding errors) have gained much attention; 
in part, this fact is related to the ever- increasing performance of computer 
hardware. When errors in initial data and rounding errors strongly affect 



Tropical Mathematics, Idcnipotcnt Analysis, Classical Mechanics, and Geometry 53 

the computation results, such as in ill-posed problems, analysis of stability 
of solutions, etc., it is often useful to perform computations with improved 
and variable accuracy. In particular, the rational arithmetic, in which the 
rounding error is specified by the user [52], can be used for this purpose. 
This arithmetic is a useful complement to the interval analysis [70]. The 
corresponding computational algorithms must be universal (in the sense 
that they must be independent of the computer representation of num- 
bers). 

13. Universal algorithms of linear algebra over semirings 

The most important linear algebra problem is to solve the system of 
linear equations 

(24) AX = B, 

where ^4 is a matrix with elements from the basic field and X and B are 
vectors (or matrices) with elements from the same field. It is required to 
find X if A and B are given. If A in (Ell) is not the identity matrix /, then 
system ( ^^ can be written in form (ESj), i.e., 

(25) X = AX^B. 

It is well known that the form (I25D is convenient for using the successive ap- 
proximation method. Applying this method with the initial approximation 
Ao = 0, we obtain the solution 

(26) A = A'B, 
where 

(27) ^* = / + ^ + ^2 ^ • • • + ^" + • • • 
On the other hand, it is clear that 

(28) A' = {I-A)-\ 

if the matrix / — ^4 is invertible. The inverse matrix (/ — A)^^ can be 
considered as a regularized sum of the formal series (1271). 
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The above considerations can be extended to a broad class of semirings. 

The closure operation for matrix semirings Mat^(5') can be defined and 
computed in terms of the closure operation for S (see Subsection 111.31 
above); some methods are described in [III21I71II1II51I2SH221I331EZII5III521IH31 
I511IH7] . One such method is described below (LDM- factorization), see |i5] . 

If 5* is a field, then, by definition, x* = {1 — x)^^ for any x ^ 1. If 5* is 
an idempotent semiring, then, by definition, 

(29) x* = lexea;^e---= sup{l, x, x^ . . . }, 

if this supremum exists. Recall that it exists if S is complete, see sec- 
tion 4.2. 

Consider a nontrivial universal algorithm applicable to matrices over 
semirings with the closure operation defined. 

Example 65 (Semiring LDM-Factorization) . Factorization of a matrix into 
the product A = LDM, where L and M are lower and upper triangular 
matrices with a unit diagonal, respectively, and D is a diagonal matrix, 
is used for solving matrix equations AX = B. We construct a similar 
decomposition for the Bellman equation X = AX © B. 

For the case AX = B, the decomposition A = LDM induces the follow- 
ing decomposition of the initial equation: 

(30) LZ = B, DY = Z, MX = Y. 
Hence, we have 

(31) A~^ = M-^D-^L-\ 

if A is invertible. In essence, it is sufficient to find the matrices L, D and 
M, since the linear system ( I5UI ) is easily solved by a combination of the 
forward substitution for Z, the trivial inversion of a diagonal matrix for Y, 
and the back substitution for X. 
Using ( 1501 ) as a pattern, we can write 

(32) Z = LZ©5, Y = DYeZ, X = MX ^Y. 
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Then 

(33) A* = MWL*. 

A triple (L, D, M) consisting of a lower triangular, diagonal, and upper 
triangular matrices is called an LDM- factorization of a matrix A if rela- 
tions ( 1521 ) and (ESj) are satisfied. We note that in this case, the principal 
diagonals of L and M are zero. 

The modification of the notion of LDM-factorization used in matrix 
analysis for the equation AX = B is constructed in analogy with a con- 
struction suggested by Carre in |Ill[I5] for Lt/-factorization. 



We stress that the algorithm described below can be applied to matrix 
computations over any semiring under the condition that the unary opera- 
tion a I—)- a* is applicable every time it is encountered in the computational 
process. Indeed, when constructing the algorithm, we use only the basic 
semiring operations of addition © and multiplication and the properties 
of associativity, commutativity of addition, and distributivity of multipli- 
cation over addition. 

If ^4 is a symmetric matrix over a semiring with a commutative multi- 
plication, the amount of computations can be halved, since M and L are 
mapped into each other under transposition. 

We begin with the case of a triangular matrix A = L (or A = M). Then, 
finding X is reduced to the forward (or back) substitution. 

Forward substitution 

We are given: 

• L = ||/*||",_i, where /* = for i < j idu lower triangular matrix with 

a zero diagonal); 
. B= ||6*||ti- 

It is required to find the solution X = ||x*||^^^ to the equation X = 
LX © B. The program fragment solving this problem is as follows. 
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foT i = 1 to n do 
{ X* := 6*; 

for j = 1 to z — 1 do 

Back substitution 
We are given 

• M = ||w!||i^o=i, where m* = for i > j (an upper triangular matrix 
with a zero diagonal); 

• B= ||6^||f=i. 

It is required to find the solution X = \\x'^\\^^i to the equation X = 
MX © B. The program fragment solving this problem is as follows. 

for z = n to 1 step —1 do 
{ x* := 6*; 

for j = n to z + 1 step — 1 do 

X* := X* (m* X*); } 

Both algorithms require (n^ — n)/2 operations and 0. 

Closure of a diagonal matrix 
We are given 

• D = diag(di,...,d„); 
. B= ||6'||ti- 

It is required to find the solution X = \\x'^\\'^^i to the equation X = 
DX B. The program fragment solving this problem is as follows. 

for i = 1 to n do 

x' ■= {d^y b'] 

This algorithm requires n operations * and n multiplications 0. 
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General case 



We are given 



• L = \\l^j\\^j^i, where /* = if z < j; 

• D = diag{di,.. .,dn)] 

• M = ||m*||". 1, where m] = if i > j: 

II J II l,J~LJ J ^ I 

. B= ||6*||ti- 

It is required to find the solution X = \\x^\\^^i to the equation X = 
AX 5, where L, D, and M form the LDM-factorization of A. The 
program fragment solving this problem is as follows. 

FORWARD SUBSTITUTION 

for z = 1 to n do 
{ X* := 6*; 

for j = 1 to z — 1 do 

x' :=x'e {l]Qx^); } 
CLOSURE OF A DIAGONAL MATRIX 
for z = 1 to n do 

x' := {d^y b'; 
BACK SUBSTITUTION 
for z = n to 1 step —1 do 
{ for j = n to z + 1 step — 1 do 

x^ := x^ (m* x-^); } 

Note that x^ is not initialized in the course of the back substitution. The 
algorithm requires n^—n operations 0, n^ operations 0, and n operations *. 

L DM- factorization 

We are given 

^ A — W^i \\n 
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It is required to find the LDM-factorization of ^4: L = ||/*||?,_i, D = 
diag((ii, . . . , dn), and M = \\m^j\\^^j=i, where /* = if z < j, and m* = if 

The program uses the following internal variables: 

• V = ||f'||^=i; 

• d. 

INITIALISATION 

for z = 1 to n do 

for j = 1 to n do 



MAIN LOOP 




fo] 


: J = 1 to 


n do 




r 
l 


for i = 


Itoj 


do 




v' 


:=a^; 





for /c = 1 to j — 1 do 

for z = /c + 1 to j do 

for z = 1 to j — 1 do 

a] := {af}* v'; 
a^- •— v^] 
for /c = 1 to j — 1 do 

for z = j + 1 to n do 

a) := a] {oJ-^Qv^); 
d = {v^Y; 
for i = j + 1 to n do 

aj := aj Q d; } 
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This algorithm requires (2n^ — 3n^ + n)/6 operations ©, {2n^ + 3n^ — 
5n)/6 operations 0, and n(n + 1)/2 operations *. After its completion, the 
matrices L, D, and M are contained, respectively, in the lower triangle, 
on the diagonal, and in the upper triangle of the matrix C. In the case 
when A is symmetric about the principal diagonal and the semiring over 
which the matrix is defined is commutative, the algorithm can be modified 
in such a way that the number of operations is reduced approximately by 
a factor of two. 



Other examples can be found in [I3[I5],[2SH2EE71EHIIH1[H7] 



Note that to compute the matrices A* and A*B it is convenient to solve 
the Bellman equation 



Some other interesting and important problems of linear algebra over 
semirings are examined, e.g., in p[I2l[I3l[ISl[23l[21[2SH22lEIl[7SH77l[72lES 



TOT 



Remark 66. It is well known that linear problems and equations are es- 
pecially convenient for parallelization, see, e.g., [97]. Standard methods 
(including the so-called block methods) constructed in the framework of 
the traditional mathematics can be extended to universal algorithms over 
semirings (the correspondence principle!). For example, formula (I22D dis- 
cussed in Subsection 11L3| leads to a simple block method for parallelization 
of the closure operations. Other standard methods of linear algebra [S^ 
can be used in a similar way. 

14. The correspondence principle for computations 

Of course, the idempotent correspondence principle is valid for algorithms 
as well as for their software and hardware implementations pOlH^lHill^T] . 
Thus: 

// we have an important and interesting numerical algorithm, then there 
is a good chance that its semiring analogs are important and interesting as 
well. 
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In particular, according to the superposition principle, analogs of linear 
algebra algorithms are especially important. Note that numerical algo- 
rithms for standard infinite-dimensional linear problems over idempotent 
semirings (i.e., for problems related to idempotent integration, integral op- 
erators and transformations, the Hamilton- Jacobi and generalized Bellman 
equations) deal with the corresponding finite-dimensional (or finite) "linear 
approximations" . Nonlinear algorithms often can be approximated by lin- 
ear ones. Thus the idempotent linear algebra is a basis for the idempotent 
numerical analysis. 

Moreover, it is well-known that linear algebra algorithms easily lend 
themselves to parallel computation; their idempotent analogs admit par- 
allelization as well. Thus we obtain a systematic way of applying parallel 
computing to optimization problems. 

Basic algorithms of linear algebra (such as inner product of two vectors, 
matrix addition and multiplication, etc.) often do not depend on concrete 
semirings, as well as on the nature of domains containing the elements of 
vectors and matrices. Algorithms to construct the closure ^4* = / © ^4 © 
A^ (B ■ ■ ■ (B A^ (B ■ ■ ■ = 0^1 A^ of an idempotent matrix A can be derived 
from standard methods for calculating (/ — ^4)^^ For the Gauss- Jordan 
elimination method (via LU-decomposition) this trick was used in [84], 
and the corresponding algorithm is universal and can be applied both to 
the Bellman equation and to computing the inverse of a real (or complex) 
matrix [I — A). Computation of A^^ can be derived from this universal 
algorithm with some obvious cosmetic transformations. 

Thus it seems reasonable to develop universal algorithms that can deal 
equally well with initial data of different domains sharing the same basic 
structure [101 HI Si] . 
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15. The correspondence principle for hardware design 

A systematic application of the correspondence principle to computer 
calculations leads to a unifying approach to software and hardware design. 

The most important and standard numerical algorithms have many hard- 
ware realizations in the form of technical devices or special processors. 
These devices often can he used as prototypes for new hardware units gen- 
erated by substitution of the usual arithmetic operations for its semiring 
analogs and by addition tools for performing neutral elements and 1 (the 
latter usually is not difficult). Of course, the case of numerical semirings 
consisting of real numbers (maybe except neutral elements) and semirings 
of numerical intervals is the most simple and natural P^liDlli2Hi^[5Tl[SHl 
[5SIED]- Note that for semifields (including Rmax and Rmin) the operation 



of division is also defined. 

Good and efficient technical ideas and decisions can be transferred from 
prototypes to new hardware units. Thus the correspondence principle gen- 
erated a regular heuristic method for hardware design. Note that to get a 
patent it is necessary to present the so-called 'invention formula', that is 
to indicate a prototype for the suggested device and the difference between 
these devices. 

Consider (as a typical example) the most popular and important algo- 
rithm of computing the scalar product of two vectors: 

(34) (x, y) = xiyi + X22/2 H h Xnyn- 

The universal version of ( jMf ) for any semiring A is obvious: 

(35) (x, y) = {xi yi) (xa 2/2) • • • {xn yn)- 
In the case A = Rmax this formula turns into the following one: 

(36) (x, y) = max{xi + yi, X2 + ?/2, • • • , a^n + Vn}- 

This calculation is standard for many optimization algorithms, so it is 
useful to construct a hardware unit for computing (BSj)- There are many 
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different devices (and patents) for computing (Elj) and every such device 
can be used as a prototype to construct a new device for computing 



and even (1351). Many processors for matrix multiplication and for other 
algorithms of linear algebra are based on computing scalar products and 
on the corresponding "elementary" devices respectively, etc. 

There are some methods to make these new devices more universal than 
their prototypes. There is a modest collection of possible operations for 
standard numerical semirings: max, min, and the usual arithmetic opera- 
tions. So, it is easy to construct programmable hardware processors with 
variable basic operations. Using modern technologies it is possible to con- 
struct cheap special-purpose multi-processor chips implementing examined 
algorithms. The so-called systolic processors are especially convenient for 
this purpose. A systolic array is a 'homogeneous' computing medium con- 
sisting of elementary processors, where the general scheme and processor 
connections are simple and regular. Every elementary processor pumps 
data in and out performing elementary operations in a such way that the 
corresponding data flow is kept up in the computing medium; there is 
an analogy with the blood circulation and this is a reason for the term 
"systolic", see e.g., HQlSZlSiiaSlESllHSlEaiHT]. 



Some systolic processors for the general algebraic path problem are pre- 
sented in [H31IB11IH7]. In particular, there is a systolic array of n{n + 1) 
elementary processors which performs computations of the Gauss-Jordan 
elimination algorithm and can solve the algebraic path problem within 
5n — 2 time steps. Of course, hardware implementations for important and 
popular basic algorithms increase the speed of data processing. 

The so-called GPGPU (General-Purpose computing on Graphics Pro- 
cessing Units) technique is another important field for applications of the 
correspondence principle. The matter is that graphic processing units 
(hidden in modern laptop and desktop computers) are potentially pow- 
erful processors for solving numerical problems. The recent tremendous 
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progress in graphical processing hardware and software resulted in new 
"open" programmable parallel computational devices (special processors), 
see, e.g., PH[751[TU2] . These devices are going to be standard for coming PC 
(personal computers) generations. Initially used for graphical processing 
only (at that time they were called GPU), today they are used for vari- 
ous fields, including audio and video processing, computer simulation, and 
encryption. But this list can be considerably enlarged following the cor- 
respondence principle: the basic operations would be used as parameters. 
Using the technique described in this paper (see also our references), stan- 
dard linear algebra algorithms can be used for solving different problems 
in different areas. In fact, the hardware supports all operations needed 
for the most important idempotent semirings: plus, times, min, max. The 
most popular linear algebra packages [ATLAS (Automatically Tuned Lin- 
ear Algebra Software), LAPACK, PLASMA (Parallel Linear Algebra for 
Scalable Multicore Architectures)] can already use GPGPU, see |103H105] . 



We propose to make these tools more powerful by using parameterized 
algorithms. 

Linear algebra over the most important numerical semirings generates 
solutions for many concrete problems in different areas, see above. 

Note that to be consistent with operations we have to redefine zero (0) 
and unit (1) elements (see above); comparison operations must be also 
redefined as it is described above. Once the operations are redefined, then 
the most of basic linear algebra algorithms, including back and forward 
substitution, Gauss elimination method, Jordan elimination method and 
others could be rewritten for new domains and data structures. Combined 
with the power of the new parallel hardware this approach could change 
PC from entertainment devices to power full instruments. 
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16. The correspondence principle for software design 

Software implementations for universal semiring algorithms are not as ef- 
ficient as hardware ones (with respect to the computation speed) but they 
are much more flexible. Program modules can deal with abstract (and 
variable) operations and data types. These operations and data types can 
be defined by the corresponding input data. In this case they can be gen- 
erated by means of additional program modules. For programs written in 
this manner it is convenient to use special techniques of the so-called object 
oriented (and functional) design, see, e.g., pUllSDll^] . Fortunately, powerful 
tools supporting the object-oriented software design have recently appeared 
including compilers for real and convenient programming languages (e.g. 
C^^ and Java) and modern computer algebra systems. 

Recently, this type of programming technique has been dubbed generic 
programming (see, e.g., P,ED])- To help automate the generic program- 
ming, the so-called Standard Template Library (STL) was developed in the 
framework of C++ fHDlEI]- However, high-level tools, such as STL, possess 
both obvious advantages and some disadvantages and must be used with 
caution. 

It seems that it is natural to obtain an implementation of the correspon- 
dence principle approach to scientific calculations in the form of a pow- 
erful software system based on a collection of universal algorithms. This 
approach ensures a working time reduction for programmers and users be- 
cause of the software unification. The arbitrary necessary accuracy and 
safety of numeric calculations can be ensured as well. 

This software system may be especially useful for designers of algorithms, 
software engineers, students and mathematicians. 

Note that there are some software systems oriented to calculations with 
idempotent semirings like Rmax! see, e.g., [^. However these systems do 
not support universal algorithms. 
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17. Interval analysis in idempotent mathematics 

Traditional interval analysis is a nontrivial and popular mathematical 
area, see, e.g., [11[211E51[7QII711|77]. An "idempotent" version of interval 
analysis (and moreover interval analysis over positive semirings) appeared 
in J5HII521ISD]- Later the idempotent interval analysis has attracted many 
experts in tropical linear algebra and applications, see, e.g., piBl [Ml 1^1751 
[TBlllOlj . We also mention the closely related interval analysis over the 



positive semiring R+ discussed in [D]. 

Let a set S be partially ordered by a relation ^. A closed interval in S 
is a subset of the form x=[x, x] = {a:G5'|x^a:^x}, where the 
elements x ^ x are called lower and upper bounds of the interval x. The 
order ^ induces a partial ordering on the set of all closed intervals in S: 
X ^ y iff X ^ y and x ^ y . 

A weak interval extension I{S) of an ordered semiring S is the set of all 
closed intervals in S endowed with operations © and defined as x © y = 
[x^y^xSy], x0y = [x0y,x0y] and a partial order induced by the 
order in S. The closure operation in I{S) is defined by x* = [x*,x*]. There 
are some other interval extensions (including the so-called strong interval 
extension [32]) but the weak extension is more convenient. 

The extension I{S) is idempotent if S is an idempotent semiring. A 
universal algorithm over S can be applied to I{S) and we shall get an 
interval version of the initial algorithm. Usually both the versions have 
the same complexity. For the discrete stationary Bellman equation and 
the corresponding optimization problems on graphs, interval analysis was 
examined in [ J5HII5D] in details. Other problems of idempotent linear algebra 
were examined in PIBIBIII^ITSIITB] . 



Idempotent mathematics appears to be remarkably simpler than its tra- 
ditional analog. For example, in traditional interval arithmetic, multipli- 
cation of intervals is not distributive with respect to addition of intervals, 
whereas in idempotent interval arithmetic this distributivity is preserved. 
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Moreover, in traditional interval analysis the set of all square interval ma- 
trices of a given order does not form even a semigroup with respect to 
matrix multiplication: this operation is not associative since distributiv- 
ity is lost in the traditional interval arithmetic. On the contrary, in the 
idempotent (and positive) case associativity is preserved. Finally, in tra- 
ditional interval analysis some problems of linear algebra, such as solution 
of a linear system of interval equations, can be very difficult (more pre- 
cisely, they are A^P-hard, see JTSKSHESIES] and references therein). It was 
noticed in [ J5S1I5D] that in the idempotent case solving an interval linear 
system requires a polynomial number of operations (similarly to the usual 
Gauss elimination algorithm). The remarkable simplicity of idempotent 
interval arithmetic is due to the following properties: the monotonicity of 
arithmetic operations and the positivity of all elements of an idempotent 
semiring. 

Interval estimates in idempotent mathematics are usually exact. In the 
traditional theory such estimates tend to be overly pessimistic. 
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